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Chapitre 1 

Introduction 



L'un des buts principaux de la Chimie est de decrire la matiere a l'echelle microscopique et correlative- 
ment de comprendre les mecanismes qui en assurent la cohesion. La rationalisation des connaissances 
en Chimie s'est longtemps heurtee a l'insuffisance des connaissances, insufhsance due aux limites des 
techniques experimentales disponibles, a l'imprecision des resultats et au succes de schemas de pensee 
qui se sont vite reveles errones. 

1.1 Necessite d'une theorie de la liaison chimique 

Pour repondre aux deux questions essentielles : 

Comment la matiere est-elle organisee ? 

Comment les corps reagissent-ils entre eux ? 
la pratique experimentale, seule, ne peut fournir qu'un catalogue de donnees qui devront ensuite etre 
classees afin de degager des tendances. Des ses origines la chimie s'est interessee a une infinie variete de 
substances complexes. Pour classer et rationaliser les connaissances, un cadre theorique est necessaire. 
Ce cadre a d'abord ete fourni par la Metaphysique, puis par la Physique. En 1699 Fontenelle ecrivait : 
<§;La Chymie par des operations visibles, resout les corps en certains principes grossiers et palpables, sels, 
soufres etc . . ., mais la Physique, par des speculations delicates, agit sur les principes comme la Chymie a 
fait sur les corps; elle les resout eux-memes en d'autres principes encore plus simples, en petits corps mus 
et figures d'une infinite de facons . . .L'esprit de la Chymie est plus confus, plus enveloppe; il ressemble 
aux mixtes, oil les principes sont plus embarrasses les uns avec les autres; l'esprit de Physique est plus 
simple, plus degage, enfin il remonte jusqu'aux premieres origines, l'autre ne va pas jusqu'au bout.^>En 
1699, bien que la theorie metaphysique des quatre elements aristoteliciens soit toujours preponderante, 
le texte de Fontenelle temoigne des rapports privilegies entre Chimie et Physique et annonce le role que 
jouera cette derniere discipline dans l'avenement d'une veritable theorie scientifique en Chimie. II suffit 
d'actualiser certains termes, de remplacer principe par element, petit corps par atome ou par electron, 
pour s'en convaincre. 

Dans ces conditions, pour que la Chimie ait le statut de science exacte, il est necessaire de disposer 
d'une theorie non-empirique de la liaison chimique. Les buts d'une telle theorie sont : 

1. de permettre une comprehension systematique des structures et des mecanismes reactionnels a 
l'echelle microscopique atomique et moleculaire 

2. de fournir des definitions rigoureuses aux termes specifiques de la chimie utilises pour decrire les 
structures moleculaires et designer les proprietes des atomes et des molecules. Afin que ces grandeurs 
soient quantifiables, il faudra etre capable d'en donner, si possible, une formulation mathematique. 

3. de permettre le calcul des proprietes des molecules et ainsi de servir de fondation a la construction 
des outils de modelisation qui sont de plus en plus utilises en recherche developpement pour creer 
des molecules ou des materiaux nouveaux ayant des proprietes specifiques. 
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L elaDoration a une tneorie satislaisant ce earner aes cnarges est line tacrie aimcile. fom s en convamcre 
il suffit de mesurer la difficulte que nous avons a definir les mots memes de liaison chimique. La definition 
la plus repandue, parce que sans doute la moins mauvaise est celle de Linus Pauling: <^.there is a chemical 
bond between two atoms or groups of atoms in case that the forces acting between them are such as to 
lead to the formation of an aggregate with sufficient stability to make it convenient for the chemist to 
consider it as an independent molecular species.^- 

Cette definition qui a les apparences d'une definition scientifique ne resiste pas a la critique. Elle 
implique tout d'abord que Ton sache definir un atome dans une molecule, de ce point de vue il aurait ete 
preferable de considerer les forces agissant sur les noyaux. II est dangereux de relier la liaison chimique 
aux forces s'exercant sur les noyaux, en effet, pour un systeme stable et done a l'equilibre, la condition 
d'equilibre impose a la valeur moyenne de la resultante des forces d'etre nulle. Si Ton s'en tient aux 
termes de la definition de Pauling il ne pourrait exister de liaison chimique au sein d'une molecule stable. 
Si les noyaux sont ecartes de leur position d'equilibre, pour la plupart des systemes polyatomiques la 
decomposition des resultantes suivant les directions internucleaires est arbitraire et il est done impossible 
de decider si deux atomes sont ou non lies. En fait dans cette definition, la decision ultime revient au 
chimiste. 

1.2 Evolution des theories 

1.2.1 La chimie aristotelicienne 

Les grecs ont ete les premiers a s'interroger sur la nature de la matiere. S'ils sont observateurs, il 
ne sont pas pour autant experimentateurs et leur reflexion reste uniquement philosophique. Pour Thales 
de Milet il existe une prima materia, l'eau, dont les differentes substances sont la manifestation de la 
diversite. Anaximandre de Milet pense que cette prima materia contient des contraires et est divisee par 
le mouvement perpetuel. Plus tard, Empedocle introduira la theorie pluraliste des quatre elements : la 
terre, l'eau, l'air et le feu. Toutes ces theories metaphysiques participent d'une description continue de la 
nature. Les partisans du discontinu, comme Leucippe et son eleve Democrite d'Abdere, introduisent le 
concept d'atome. Les atomes sont des solides homogenes, en nombre d'espece fini et en nombre fini. lis 
sont spontanement animes de mouvement et leurs interactions quand ils se deplacent dans le vide forment 
le monde sensible. 

Les eleves de Socrate, Platon et Aristote ne partagent pas le point de vue atomiste. Pour eux, il existe 
cinq elements, dont quatre terrestres, auxquels correspondent les cinq polyedres reguliers. II existe de plus 
pour Aristote, quatre qualites : le chaud, le froid, le sec et l'humide. Chaque element terrestre possede deux 
qualites. La transmutation d'un element en un autre est possible par echange de deux qualites. La theorie 
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Tab. 1.1: Correspond an ce entre elements, polyedres reguliers et qualites. 

des quatre elements va durer sous differentes formes jusqu'a la fin du XVIII 6 siecle. A la Renaissance, 
Paracelse enonce une theorie ou les quatre elements sont completes par trois principes : le sel, le soufre 
et le mercure. Ces trois principes ont pour origine l'alchimie arabe et decrivent le comportement de la 



matiere vis a vis cm leu : le souire est associe a la comDustion, le mercure ail cliangement a etat pliysique 
et le sel a l'inertie. 

1.2.2 La chimie newtonienne 

La theorie des quatre elements est attaquee au debut du XVII 6 siecle, d'abord par van Helmont qui 
le premier invente le concept de gaz apres avoir prepare du dioxide de carbone (gaz sylvestre) et des 
melanges d'hydrocarbures (gaz pingue). II distingue ces gaz de l'air ordinaire montrant ainsi l'insuffisance 
de cet element generique. 

Rene Descartes decrit la matiere comme un ensemble de corpuscules en mouvement susceptibles de 
se transformer, de se regrouper, de se diviser sous l'effet des chocs. II n'existe pas de force d'interaction 
entre ces particules en raison de l'absence d'un support a une telle force. La constitution de produits 
a partir de reactifs est expliquee en invoquant des facteurs de forme, on invoque des atomes crochus 
pour expliquer la formation des mixtes, ainsi pour Lhemery les acides sont des solides munis de pointes 
pouvant s'emboiter dans les trous des bases. 

La premiere theorie explicative de la liaison chimique est due a Isaac Newton. Newton s'est surtout 
interesse a la chimie et a rechercher un principe expliquant pourquoi les substances chimiques reagissent 
entre elles, cela l'a amene a etudier le mouvement des planetes et a formuler la theorie de la gravitation 
universelle. Avant Newton, le mouvement des planetes est explique en invoquant une force d'interaction 
entre le soleil et chaque planete. Pour Newton il faut considerer le systeme solaire dans son ensemble et 
tenir compte de l'interaction mutuelle des astres qui s'exprime a l'aide d'une force universelle en 1/r 2 
et qui est proportionnelle aux masses des astres. Newton pense que les lois de la gravitation universelle 
s'appliquent egalement a la Chimie, il faut considerer les systemes dans leur ensemble, ainsi quand un 
solvant rend une reaction possible, des particules insociables ont ete rendues sociables par la mediation 
d'un tiers. 

La force qui agit entre les corps chimiques est l'affinite, ainsi definie par Berthollet :<§;Les puissances 
qui produisent les phenomenes chimiques sont toutes derivees de l'attraction mutuelle des molecules a 
laquelle on a donne nom d'affinite, pour la distinguer de l'attraction astronomique. II est probable que 
l'une et l'autre ne soit qu'une meme propriete.^> 

Le modele newtonien est perfectionne par Roger Boscovitch qui en 1758 publie la <§;Theorie de la 
Philosophic Naturelle^>ou il assimile les corps chimiques a des masses ponctuelles formant des edifices 
compliques. II pense que la loi d'attraction en 1/r 2 n'est valable qu'aux grandes distances. Aux courtes 
distances les forces peuvent etre attractives ou repulsives. Les forces d'attraction qui caracterisent les 
rapports entre deux corps sont specifiques parce que determinees par 1 'edifice, elles sont la resultante des 
forces dont chaque point est porteur. 

Pour Buffon les corps chimiques ne peuvent etre reduits a un systeme de masses ponctuelles et la force 
d'attraction etant en 1/r 2 , la specificite et la diversite des actions chimiques s'expliquent par la diversite 
des figures des corps. II pense qu'un jour on sera capable de deduire les affinites a partir des figures des 
particules et de prevoir les reactions par le calcul. 

Pour decouvrir les rapports entre les substances, les chimistes du XVIII 6 siecle etudient les reactions 
de deplacement dans lesquelles un acide ou une base est remplace par un autre acide ou une autre base 
dans un sel. Les premieres tables d'affinite ont ete etablies par Geoffrey en 1718. T. Bergman, cinquante 
ans plus tard fait entrer les conditions de reaction dans les tables d'affinite et introduit le concept d'affinite 
elective. 

La theorie de l'affinite ne s'est pas substitute a la theorie des quatre elements qui sert de base a 
la theorie du phlogiston developpee par Stahl pour expliquer la combustion. Cette theorie forme un 
systeme coherent qui est capable d 'expliquer qualitativement la plupart des connaissances chimiques du 
XVIII 6 siecle. Cette theorie conserve les quatre elements aristoteliciens et introduit un principe sans 
masse, le phlogiston qui est perdu par les corps au cours des combustions et qui est acquis par reduction. 
Pour Kant, Stahl est un heros de la methode experimentale. En 1766, Pierre Joseph Macquer ecrit dans 
le <CDictionnaire de Chymie^> : <§;On reconnaitra sans doute avec etonnement que nous admettons a 
present comme principes de tous les composes les quatre elements, le feu, Fair, l'eau et la terre, qu'Aristote 



avait maiques comme tels avant qu on eut les connaissances necessaires pour constater pareille vente.;j> 

1.2.3 Dalton et le renouveau de l'atomisme 

Le coup fatal a la theorie des quatre elements est porte a la fin du XVIII 6 siecle par Lavoisier. Paral- 
lelement a la mise en place du systeme metrique, un effort de rationalisation est entrepris en Chimie par 
l'elaboration d'une nouvelle nomenclature. En etablissant cette nomenclature, Lavoisier et ses collegues 
Fourcroy, Guyton de Morveau et Berthollet rejettent les quatre elements classiques et publient une table 
comprenant trente trois elements, dont la lumiere et le calorique. Les elements entrent dans la composition 
des corps purs, un seul element participe a un corps simple. 

Dalton adhere a cette nouvelle nomenclature, et revitalise la theorie atomiste. A chaque element 
correspond un atome constitue d'un petit corps solide entoure d'une enveloppe de calorique. En 1803, 
il publie une premiere table de masses atomiques. La theorie atomiste de Dalton trouve des justifica- 
tions experimentales dans la loi des proportions definies. Cette theorie va permettre de proposer des 
stcechiometries pour les differents composes. 

Comme l'a note G. N. Lewis, la loi des proportions definies a transforme une speculation philosophique 
en une theorie scientifique feconde. La theorie atomiste et moleculaire est a l'origine non seulement de la 
chimie structurale mais aussi de la thermodynamique statistique. Les travaux de Dalton ont fait evolue 
la conception de la matiere d'une vision continue a une description discontinue. La fin du XIX e siecle 
apportera les preuves experimentales directes de la discontinuite de la matiere. 

La theorie atomiste de Dalton rend compte de la loi des proportions multiples, mais ne permet pas 
de comprendre la loi de Gay-Lussac de combinaison des gaz en volume. Si Ton suppose comme Berzelius 
que des volumes egaux de gaz contiennent le meme nombre d'atomes (atome signifiant la plus petite 
particule d'un compose) on aboutit a une absurdite pour un grand nombre de reactions. Dans le cas de 
la combinaison de l'hydrogene et du chlore 

H + CI -► HC1 

l'hypothese de Berzelius suggere que les atomes d'hydrogene et de chlore se divisent en deux ce qui est 
contraire au point de vue de la theorie atomiste. 

Cette difficulte a ete surmontee en 1811 par Amedeo Avogadro qui introduit la notion de molecule. 
Une molecule est la plus petite quantite de matiere physiquement isolable, c'est a dire existant a l'etat 
libre en phase gazeuse. Les molecules sont constitutes d'atomes isoles ou de groupements d'atomes lies. 

La conjecture de Prout est une autre idee feconde du XIX e siecle. En 1815, en examinant les valeurs 
les masses atomiques deja proposees Prout remarque qu'elles sont a peu pres des multiples en tiers de la 
masse atomique de l'hydrogene. II en deduit que l'hydrogene est la primera materia a partir de laquelle 
les differents elements sont constitues. 

Avec la connaissance de plus en plus precise d'un plus grand nombre de masses atomiques des ecarts 
significatifs a la loi de Prout ont ete mis en evidence. Cependant, un calcul de probability montre que 
la loi de Prout ne peut etre accidentelle. Rydberg a calcule en 1887, qu'il n'y a qu'une chance sur un 
milliard pour que les vingt deux premiers elements aient des masses atomiques aussi proches de nombres 
entiers. 

1.2.4 L'evolution des idees au XIX e 

Au XIX e siecle on assiste a un essor considerable de la Chimie. Du point de vue de la Chimie generale 
le progres le plus significatif est realise par Berzelius qui introduit le systeme de notation moderne et qui 
determine avec une grande precision les masses des equivalents qui deviendront par la suite les masses 
atomiques. 

Les experiences de Davy (1778-1829) sur les effets du courant l'electrique sur les substances chimiques 
suggerent que l'electricite joue un role primordial dans la liaison chimique. Pour lui l'affinite chimique est 
de nature electrostatique. A partir de cette idee, Berzelius developpe en 1819 une theorie electrochimique 
(ou dualistique) de la matiere qui sera la theorie dominante de la premiere partie du XIX e siecle. D'apres 
la theorie dualistique, lorsqu'un atome entre en contact avec un autre atome plus electronegatif, il se 



produit im nux a electncite tel que 1 im aes atomes se cliarge positivement et 1 autre negativement. La 
cohesion de l'edifice moleculaire ainsi forme est alors due aux forces electrostatiques attractives entre 
atomes portant des charges opposees. Comme on peut le remarquer, cette theorie est tres proche du 
modele contemporain de la liaison ionique. 

La theorie dualistique est appliquee avec succes aux sels et aux composes mineraux. Quand la theorie 
electrochimique fut proposee, on ignorait encore que beaucoup de molecules fortement liees sont consti- 
tutes d'atomes identiques comme c'est le cas du dihydrogene et du diazote. L'existence de telles molecules 
est en contradiction avec les hypotheses de la theorie dualistique qui suppose une asymetrie de la liaison. 
Le developpement de la chimie organique a egalement contribue a la critique de la theorie dualistique. En 
particulier, le fait que l'hydrogene repute electropositif puisse etre substitue par un halogene electronega- 
tif sans changer notablement les proprietes d'un compose organique, ne peut etre explique par la theorie 
electrochimique . 

Faraday a montre en 1833 que la loi des proportions multiples s'applique non seulement aux masses et 
aux volumes mais egalement a la charge electrique. II a fallu plusieurs decennies avant que Ton ne realise 
que la loi de Faraday implique la discontinuite de la charge electrique. Helmholtz, le premier en 1881, 
emit l'hypothese de l'existence d'un <Catome^> d'electricite qui sera confirmee experimentalement par la 
decouverte en 1897 de 1 'electron par Joseph John Thomson. Les experiences de Thomson ont permis de 
mesurer le rapport e/m. La charge de l'electron sera mesure par la suite par Thomson, H. A. Wilson et 
Millikan. II est interessant de constater qu'il etait possible de calculer la charge de l'electron des 1865 a 
partir de la loi de Faraday et du nombre de Lodschmidt. 1 

Le developpement de la Chimie organique a d'autres consequences tres importantes : l'introduction 
du concept de valence par Frankland et Kolbe vers 1850, puis la mise au point de formules developpees 
par Couper et Kekule. Plus tard, la decouverte d'isomeres optiques par Pasteur conduira Le Bel et Van't 
Hoff a formuler l'hypothese du carbone tetraedrique : la structure moleculaire moderne est nee. La theorie 
de la stereochimie est etendue aux complexes inorganiques par Werner. 

Les reflexions sur la valence, qui est supposee alors avoir une valeur unique et caracteristique pour 
chaque element, sur les masses atomiques, sur les similarites de proprietes physico-chimiques, menees 
dans le contexte de l'explosion demographique des elements du XIX e siecle se concretiseront dans la 
classification periodique des elements de Mendeleev. 

1.2.5 Les modeles prequantiques de l'atome au debut du XX e siecle 

Apres la mise en evidence de l'electron par Joseph John Thomson en 1897, un grand nombre de modeles 
de l'atome ont ete proposes. Parmi ceux-ci nous decrirons le modele de Thomson lui-meme datant de 
1904, ou les electrons sont immerges a l'interieur d'une sphere uniformement chargee d'electricite positive 
et celui de Lewis ou ils sont disposes aux sommets d'un cube. Dans le modele de Thomson les electrons 
gravitent sur des anneaux, lorsque le nombre d'electrons d'un anneau devient superieur a une limite 
donnee, cet anneau se divise en deux anneaux concentriques. Ce modele rend ainsi compte des proprietes 
atomiques issues de la classification periodique. Dans son modele de 1902, non publie, Lewis formule 
quatre hypotheses principales: 

1. Les electrons sont disposes suivant des cubes concentriques 

2. L'atome neutre d'un element possede un electron de plus que l'atome neutre de l'element precedent 

3. Le cube a huit electrons est atteint pour les atomes de gaz rare, ce cube est le noyau autour duquel 
s'organise le plus grand cube electronique de la ligne suivante 

J En 1865, Loschmidt calcule le nombre de molecules contenues dans 1 cm 3 de gaz a 0° C et sous 1 atmosphere a partir 
de la theorie cinetique des gaz. 

£ = 2.691 10 19 molecules cm -3 
e = TjC X 22414 = 96500/2.691 10 19 X 224141 
= 1.5999 10 _19 C 

Cette valeur est a comparer a la valeur des tables: 1.6021917 10 — 19 C. 



4. les electrons a im cuDe externe mcomplet peuvent etre cedes a im autre atome arm qu 11 complete 
son cube externe a huit electrons. 

La premiere publication faisant etat de la stabilite du groupe de huit electrons est due a Abegg en 1904. 

En 1911 Rutherford decouvre l'existence du noyau. L'atome devient une structure ouverte constitute 
d'un noyau charge positivement et d 'electrons. Nagaoka a montre, des 1904, qu'une telle structure est 
stable si l'attraction electron-noyau est forte et si les electrons sont animes d'un mouvement rapide. 

En 1916, G. N. Lewis ameliore son modele de couches cubiques maintenant concentriques autour du 
noyau. La premiere couche contient deux electrons. La meme annee, Kossel propose un modele equivalent 
constitue d'orbites circulaires concentriques. 




Lewis Kossel 

Fig. 1.1: Representation de l'atome d'argon suivant Lewis et Kossel 

En 1913, Bohr propose le premier modele quantique pour l'hydrogene qui donne une interpretation 
qualitative du spectre atomique. 



1.2.6 Theorie de Lewis de la valence 

La connaissance de la structure de l'atome montre que seules les forces electrostatiques peuvent etre 
invoquees pour expliquer la stabilite des molecules. De ce point de vue la theorie dualistique de Berzelius 
etait satisfaisante, mais il s'agit maintenant d'en corriger les insuffisances. En 1907, J.J. Thomson a 
suggere que deux atomes neutres pouvaient etre lies par des forces electrostatiques sans qu'il y ait transfert 
d'electricite de l'un vers l'autre. Pour cela, il considere deux spheres identiques d 'electrisation uniforme 
se recouvrant partiellement. En 1908, Ramsay pense que des electrons peuvent etre partages entre deux 
atomes. Cette approche est developpee par Stark dans le cadre de sa theorie dynamique atomique. Dans 
ce modele, les electrons sont localises entre les deux atomes et forment la liaison chimique. 

L'idee de Lewis est que dans une molecule les atomes tendent a saturer leur couche de valence a huit 
electrons, non seulement par transfert d'un atome a l'autre, mais aussi en partageant une ou plusieurs 
paires d'electrons. Ces electrons partages sont consideres comme appartenant aux couches externes de 
chacun des atomes. L'originalite de Lewis a ete de substituer une regie de paire a la regie de 1 'octet ce qui 
permet de rendre compte de la structure electronique de l'helium et des composes hydrogenes. D'autre 
part, la plupart des molecules stables ont un nombre pair d'electrons. La representation que Lewis se fait 
de la liaison chimique est celle de deux electrons apparies appartenant aux couches externes des deux 
atomes lies. Dans les schemas de Lewis, chaque electron de valence est represente par un point: 
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h ; h 




h ; c ; h 



H 



La polarite des liaisons est representee en deplacant le doublet electronique d'un centre atomique vers 
l'autre. 

Cette intuition geniale des paires d'electrons trouvera sa justification plus tard avec la decouverte du 
spin electronique et des proprietes des fermions. La theorie de Lewis explique qualitativement l'existence 
des liaisons multiples. Dans cette theorie, la valence d'un atome dans une molecule est egale au nombre 
de paires qu'il partage avec les autres atomes. La theorie de Lewis rend compte des exceptions a la regie 
de l'octet tant pour les atomes hypervalents comme S dans SF$ que pour les atomes trivalents comme le 
bore. 

Cette theorie a apporte de nouvelles idees non seulement en ce qui concerne la structure moleculaire 
mais aussi sur la reactivite comme l'atteste la definition des bases et des acides proposee par Lewis. Une 
base est une substance qui possede une paire d'electrons susceptible d'etre cedee a un autre atome afin 
d'en completer la couche externe. Un acide est une substance susceptible d'utiliser une paire d'electrons. 

1.2.7 Importance de la Mecanique quantique 

Le modele quantique de l'atome propose en 1913 est le premier a donner une base non-empirique 
a la formule de Balmer. II permet de prevoir avec une tres grande precision les longueurs d'onde des 
transitions electroniques de l'atome d'hydrogene et des ions hydrogenoi'des, par la il s'avere predictif. 
La serie de Pickering, observee dans le spectre de l'etoile (*-Puppis qui avait ete attribute par erreur a 
l'hydrogene est attribute a He+ . Le modele de Bohr est ameliore en 1918 par Sommerfeld qui introduit 
des orbites elliptique afin de rendre compte de la structure fine des niveaux d'energie de l'hydrogene et 
des ions hydrogenoi'des. Les model es de Bohr et de Sommerfeld sont limites a des cas particulier bien 
qu 'important et de ce fait ils ne sauraient servir de base a une description quantique de la liaison chimique. 

Le fondement axiomatique de la mecanique quantique est apporte par Louis de Broglie en 1924 avec 



qui postule la dualite onde corpuscule. Cette relation est le point d'appui de l'equation de Schrodinger 
publiee en 1926. Avec l'equation de Schrodinger completee par le principe de Pauli il devient, en principe, 
possible de traiter exactement de facon non relativiste tout systeme microscopique de particules. L'equa- 
tion de Dirac qui est l'equivalent relativiste de l'equation de Schrodinger introduit naturellement, c'est a 
dire sans en postuler l'existence, le spin de l'electron. Les premieres applications aux molecules datent de 
1926 en ce qui concerne l'ion moleculaire traite par Burrau et de 1927 pour la molecule de dihydrogene 
etudiee independemment par Condon et par Heitler et London. Dans le cas du dihydrogene, il n'existe 
pas de solutions exactes a l'equation de Schrodinger ce qui contraint a faire des approximations. L'ap- 
proximation des orbitales moleculaires utilisee par Condon et l'approximation valence-bond d'Heitler et 
London sont devenus les points de depart methodologiques de toute etude de systemes polyelectroniques. 
Elles sont aussi a l'origine des theories qualitatives et explicatives de la liaison chimique: theorie de la 
formation des orbitales moleculaires a partir d'orbitales atomiques, theories de Pauling de l'hybridation 
et de la resonance. Des solutions approchees de bien meilleur qualite ont ete proposees, d'abord pour 
l'helium par Hylleraas puis en 1933 pour le dihydrogene par James et Coolidge. Bien que de tels calculs 
permettent d'atteindre la precision spectroscopique, ils n'ont pas suscite de generalisation qualitative. 

Du point de vue de la Chimie quantique quantitative, les progres de l'electronique et de l'informatique 
ont entraine la mise au point de puissants logiciels de calcul des structures electroniques. 



la relation 



A = h/mv 



Un progres tneonque tres important a eii lieu en iyo4 avec le tneoreme de rionenDerg et Kolin qui 
est le theoreme d'existence de la fonctionnelle de la densite. La theorie de la fonctionnelle de la densite 
fournit une alternative a la resolution de l'equation de Schrodinger. 

Les theories explicatives, evoquees precedemment souffrent d'apres moi, du fait qu'elles derivent d'une 
description avec les mots de la Physique de modes operatoires, eux-memes tributaires d'approximations. 
Paradoxalement, elles sont mal ad ap tees a la description de la fonction d'onde exacte ou meme a celle 
d'une fonction d'onde approchee de bonne qualite. La theorie des loges de R. Daudel est beaucoup plus 
rigoureuse, mais difficilement quantifiable, elle n'a pas eu le rayonnement qu'elle merite. La theorie de 
Bader des atomes dans les molecules est une theorie explicative rigoureuse qui reconcilie le point de vue 
localise de la chimie et la delocalistion propre a la mecanique ondulatoire. 

Le vieux projet de Newton et de Buffon se realise grace a la mecanique quantique. On dispose d'une 
theorie physique qui permet en principe d'expliquer toute la chimie microscopique et de faire des calculs 
predictifs. 



Deuxieme partie 

Mecanique Quantique des systemes 

polyelectroniques 
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Chapitre 2 



L'equation de Shrodinger et la 
fonction d'onde. 



L'application de la mecanique quantique aux systemes constitues par plusieurs particules comme 
les atomesi, les molecules ou les cristaux sera aborde dans ce chapitre dans le cadre non-relativiste de 
l'equation de Schrodinger. Comme nous le verrons, il s'agit principalement d'une generalisation du cas a 
une particule completee par un ensemble de conditions rendant compte de la nature physique du systeme 
considere. 

2.1 Construction de l'operateur hamiltonien 

Une molecule est un systeme quantique constitue de M noyaux de masse de charge +Zj^e reperes 
par les vecteurs position et de N electrons de masse m e , de charge — e et de position r,. Dans le cas 
le plus general, l'hamiltonien electronique d'un systrhe moleculaire s'ecrit sous la forme d'une somme de 
contributions : 



Dans cette expression Tjy et T e sont les operateurs energie cinetique des noyaux et des electrons qui 
s'expriment tout deux sous la forme de sommes de contributions individuelles : 



Les operateurs energie potentielle electron-noyau, noyau-noyau et electron-electron sont les sommes des 
interactions coulombiennes entre particules chargees : 



H = T N + T e + V e>e + V e>N + V N>N + V N>ext + V e , 



(2.1) 




(2.2) 




(2.3) 




(2.4) 




(2.5) 



e,e 




(2.6) 
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linnn Vjv.ext et Ve,ext sont les operateurs ae 1 energie potentielle a interaction aes noyaux et des electrons 
avec le champ electromagnetique externe dont nous ne donnons ici que l'expression formelle : 



M 



V N , ext = ^ V ext (R A , a A ) (2.7) 

A 
N 

V e ,ext = Y,^^^^) ( 2 - 8 ) 



afin de faire apparaitre les coordonnees de spin a A et 07 ainsi que le caractere additif de ces potentiels. 
Meme si Ton neglige la plupart du temps les potentiels de couplage avec un champ externe, leur expression 
montre que la fonction d'onde totale solution de l'equation de Schrodinger est une fonction de 4(M +N) + l 
variables : 3(M + N) + 1 variables d'espace, (M + N) variables de spin et le temps : 

^{R^^K^},*) (2.9) 

II est commode de regrouper les coordonnees d'espace et de spin de chaque particule dans un vecteur a 
quatre composantes qui sera designe par s'il s'agit d'un noyau et par Xj s'il s'agit d'un electron. 

R A \ ( r,: 



X - = U"J< X ' = UJ (210) 

L'equation de Schrodinger dependant du temps a alors pour expression : 

H^({X A },{^},t) = i h^({X A },{ Xi },t) (2.11) 

Toutes les fonctions Psi({X^}, {x,}, t) solution de l'equation de Schrodinger dependant du temps ne sont 
pas acceptables du point de vue de la physique. L 'interpretation statistique de la mecanique quantique 
implique que \f*({X J i}, {x,}, ^^({X^}, {xj},i) represente la probabilite de trouver le systeme au temps 
t dans la configuration de la donnee des coordonnees {X^j^jx;}, en d'autre termes la fonction d'onde 
doit avoir une valeur finie en tout point de l'espace a 4(M + N) dimensions. L'origine des principales 
autres conditions sont la nature des particules, le spin et la symetrie du potentiel nucleaire. 



2.2 L'indiscernabilite des electrons et le principe d'antisymetrie 

a quatre composantes note x,: La fonction d'onde electronique du systeme est une fonction des va- 
riables X,: 

*(X1,X 2 , ...,Xi,...,Xj,. . .,xjv) 

. Comme les electrons sont indiscernables, la permutation de deux quelconques d'entre eux ne modifie pas 
la situation physique du systeme. Cette permutation est realisee en appliquant l'operateur de transposition 
Pij a la fonction d'onde electronique. L'operateur de transposition permute avec l'operateur hamiltonien 
parceque d'une part l'addition est commutative et que d'autre part |rj — r j | = |rj — rj\. Les fonctions 
propres de H doivent etre egalement fonctions propres de Py. Designons par a la valeur propre de P,j 
associee a la fonction propre \f(xi , X2, . . . , Xj, . . . , Xj, . . . , xjv). 

P ij *(xi,X 2 , . . .,Xj, . ..,x j: . . .,Xjv) = *(xi,x 2 , ...,Xi,. • • .,x,-, . . .,Xjv) 

= a*(xi,x 2 ,...,x i ,...,Xj,...,x JV ) (2.12) 

Les transposition P,j sont des elements d'un groupe, appele groupe des permutations et note Sjv, plus 
precisement ce sont les generateurs du groupe. Comme Pjj est son propre inverse P;j x Pij — E, oil E 
est l'identite, la valeur de a peut etre determinee en appliquant deux fois Py-: 

P i P i *(xi,x 2 , . . .,xi, . ..,Xj, . . .,xjv) = a 2 M?(xi,X2, . . .,x,, . ..,Xj, . ..,x N ) 

= *(xi,x 2 , . . .,xi, ■ ■ .,Xj, . . .,xjv) (2-13) 



qui implique: 

a = ±l (2.14) 

La valeur propre +1 correspond aux systeme de spin entier, les bosons, tandis que la valeur propre —1 
correspond aux fermions (particules de spin demi-entier). Les electrons sont des fermions et la fonction 
d'onde electronique est done antisymetrique par rapport a la permutation des coordonnees d'espace et 
de spin de deux electrons quelconques: 

Pij*(xi,x 2 , . . .,xi, . ..,Xj, . . .,xjv) = -*(xi,x 2 , . . .,xi, . ..,Xj, . . .,xjv) (2-15) 

L'equation 2.15 exprime le principe d'antisymetrie ou principe de Pauli. Une consequence importante est 
que deux electrons de spins paralelles ne peuvent pas occuper la meme position. Supposons Xj = x,, alors 

Pij*(xi,X 2 , ...,Xi,. ..,Xj, . . .,Xjv) = *(X1,X 2 , ...,Xi,. ..,Xj, . ..,Xi,X N ) 

= -*(xi,x 2 , . . . ,Xj, . . .,x,,xjv) (2-16) 

relation qui ne peut etre verifiee que si: 

*(xi,x 2 , ...,Xi,...,x j: . . .,xjv) = (2.17) 
Par contre pour des electrons de spins antiparalelles 

%)> x i=( ) ' donc ' x i * Xi 

la fonction d'onde ne s'annule pas obligatoirement parceque la divergence intoduite par l'operateur V eie 
est exactement compensee par l'energie cinetique. Physiquement tout se passe comme si les electrons 
de spins paralelles etaient soumis aux coutes distances a un potentiel repulsif plus fort que le potentiel 
coulombien classique. Cette contribution purement quantique est appelee repulsion de Pauli. 

2.3 Relations d'Heisenberg et d'Ehrenfest 

La derivee par rapport au temps de la valeur moyenne de l'observable A est donnee par la relation 
d'Heisenberg : 

d{A)/dt=(%/h){W\[H,A]\W) + {dA/dt) (2.18) 
En derivant (^1^41^) par rapport au temps on obtient : 

d 3* - - 3* dA 

^<*W*> = (-^1^1*) + W A \-Q[) + <*l -QfW ( 2 - 19 ) 

Les derivees du bra et du ket sont evaluees directement a partir de l'equation de Schrodinger dependant 
du temps 

H\y) = ih\ — } (2.20) 



et 



at 



d i i f) A 

-(A) = -(*|^|*>--(*|^|*> + (*| — 1*> (2.21) 
= L(y\[H,A]\V) + (V\^\V) (2.22) 



Une observable B telle que 



[H,B] = 0, 8B/dt = (2.23) 



est appelee constante cm mouvement. 

La premiere relation d'Ehrenfest est obtenue en choisissant A = r : 



m[H,r] = -(-thV) = -p (2.24) 

i i 

et 

m|(r> = (p) (2.25) 

Avec ^4 = p 

[If, p] = z/i W(r) = -iftF(r) (2.26) 

et 

|(P) = <F(r)> (2.27) 

dans cette expression, la deuxieme relation d'Ehrenfest, F(r) est l'operateur correspondant a la force. 
En combinant ces deux relations, on obtient une equation analogue a l'equation de Newton : 

= (F(r)) (2.28) 

2.4 Theoreme du viriel 

En appliquant l'equation d'Heisenberg a l'operateur 

i = r-p (2.29) 

on obtient : 

^<r-p) = -l(|[jy,r-p]|) (2.30) 

avec 

i - 7i 2 - 

-[Jf, r • p] = ( V 2 - r • W) = 2T - r • VV (2.31) 

n m 

ou T designe l'operateur energie cinetique. 

^(r-p) = 2(f> + (-r-VV> (2.32) 

pour un etat stationnaire ou 

^■P> = (2.33) 

on obtient le theoreme du viriel 

2(f) - (r • VV) = (2.34) 

Si l'operateur energie potentielle V est une fonction homogene de degre n des coordonnees {r} le theoreme 
d'Euler 1 permet d'ecrire 

r-VV = nV (2.35) 

et 

2(f) = n(V) (2.36) 
Pour un atome isole l'operateur V est de degre -1 

2(f) = -(V) (2.37) 



1 F(x, y, z) est une fonction homogene des variables x, y, z de degre n si 

f(tx, ty, tz) = t n F{x, y, z) 

Alors 

dF dF dF . 

x ~. Hi/— V z — = nF(x,y, z) 

ox oy oz 
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Si l'hamiltonien d'un systeme depend d'un parametre A alors : 

A £ ( A ) = <*(A)|^H|*(A)) (2.38) 

L'equation de Schrodinger implique : 

(¥(A)|#(A) - £(A)|¥(A)) = (2.39) 

done 

J ; <*(A)|jj(A)-.E(A)|*(A)) = (^ffil|A(A)-£(A)|*(A))+(*(A)|a(A)-£(A)|^l) + 
comme 

(^(A)- J B(A))|*(A)>=0 (2.41) 
que E(X) est un scalaire et en choisissant |\t(A)) normee 

±E{\) = (xp(A)|^^|xp(A)} (2.42) 

2.6 Theorie des perturbations 

La methode de perturbation permet de calculer les fonctions propres et les valeurs propres approchees 
d'un hamiltonien H quelconque mais voisin d'un hamiltonien if (°) dont on connait les fonctions propres 
et les valeurs propres. En designant par A un parametre arbitraire, appele parametre de perturbation et 
qui sera suppose petit, if s'ecrit : 

H = H (0) + XH (1) (2.43) 

if et etant hermitiques, l'est egalement. Si if' 1 ) est independent du temps, la theorie des per- 
turbation est dite statique. Si depend du temps, elle est appelee theorie des perturbations dependant 
du temps ou perturbation dynamique. 

2.6.1 Perturbations statiques. 
Etats non degeneres 

Pour resoudre l'equation de Schrodinger dans le cas d'un etat non degenere 

(H - Ei)\Vi) = {HW + XH^ - Ei)\Vi) = (2.44) 
on developpe E{ et en serie de A 



Ei = E\ 0) + XEP + X 2 E^ + ... = £ XkE i k) ( 2 - 45 ) 

k = 



|*i> = |*| 0) >+A|*j 1) > + A 2 |*j 2) > + ... = ^A*|*j* ) > (2.46) 

k = 

en substituant dans l'equation 2.44, on obtient 

oo / k \ 

J2X k HW\VW) + H^\^j k - 1) ) - = (2.47) 
k=i \ j=o I 



pour que cette equation soit vermee maepenaamment de la valeur de A, 11 taut que tous les coemcients 
de X k soit nuls, ce qui engendre la suite d'equations 



(£(o) _ + (#U) _ ijW)!^- 1 )) -^EP\^ k - j) ) = (2.48) 



j = 2 



La normalisation de |^,-°^} et de 1^,} implique 2 



oo k 

j(j)|lT/(' ! -j)\ 



x;a*x)<*.- i*.- >=° ( 2 - 49 ) 

* = 1 i = o 



en multipliant a gauche l'equation 2.48 par (M*> 



(°)| 



(^ItfU)^- 1 )} = ^£iW)(*( )|*(*-j)> (2.50) 
i=i 



comme est hermitique 



=^^)(*(*-^|*(°)> (2.51) 
i=i 



done 



<^ 0) |^ (1) l^f- 1} > = ^^^^{(^f^^-^) + {^r-^l^ ^} (2-52) 

en particulier 



2 

i=i 



EP = (vfWW) (2-53) 
£< 2) = (^ 0) \H^W ] ) (2-54) 

La resolution analytique des equations 2.48 n'est pas toujours possible. Pour resoudre on utilisera le 
fait que les fonctions propres de l'operateur hermitique forment une base complete sur laquelle on 

peut projeter En introduisant le projecteur ^ |\P; C ^)(\P; C ^ | dans 2.48 il vient 



£ (£<°> - 4 0) )l< ) }(^ (0) l^ ) ) + (1) - £; i (1) )|* 1 (0) >(* 1 (0) |* i ( *- 1) > - E^I^X^N*^} 

l [ j = 2 

(2.55) 

Au premier ordre 

- ^)|*W)(*W|*^> + l^}^^ 1 )^} = (2.56) 



2 Un autre type de normalisation est tres utilise. II impose que les corrections a la fonction d'onde ) soient orthogo- 

nales a |*( C) ). 

)*(*)) = 

dans ces conditions 

E W = ^("JiHt 1 )^*- 1 ') 



done 



et 



E^ = 



,(°)\ 



(0) P (0) 



(2.57) 



(2.58) 



Etats degeneres 

Pour un etat d fois degenere, les denominateurs des equations 2.57 et 2.58 peuvent etre nuls sans 
que les numerateurs le soient. Dans ce cas la correction de la fonction d'onde au premier ordre ainsi que 
la correction du deuxieme ordre a l'energie divergent, et ne sont done pas physiques. Deux methodes 
peuvent etre utilisees pour traiter le cas des niveaux degeneres. 

Dans la premiere methode, on recherche une transformation unitaire U des fonctions I'S^}, avec 
a= l,...,d 

$ = #U 



telle que la matrice 



HW($) = 



/ (*ii|^ (1) |**i> ••• {$ii\H^\$ id ) \ 



V {^H^a) ... {<S> id \H^\<S> id ) J 
soit diagonale. U est la matrice des vecteurs propres de la matrice 

/ OM^Ml^i) ... <¥ii|ff (1) |¥tt) \ 
HW(tf) : : : 

V 0M# (1) l^i) ... (V id \HW\V id ) J 

En utilisant les fonctions on obtient 



l*&> 



(2.59) 



(2.60) 



(2.61) 



(2.62) 
(2.63) 



et 



Les termes correspondant a e\°^ = E^ sont exclus des sommes. 
-ffW est diagonal sur la base des done 



(2.64) 



8=1 



(2.65) 



En substituant dans l'equation 2.62, en multipliant par les coefficients Up a , en sommant sur l'indice /? et 
en utilisant l'egalite des projecteurs 



(2.66) 



11 vient 

(*<°> - e^) e "iM?) + E «*«i*, ( ?> - E E «ki*, ( S Wi* (1 K } > = 

/3 = 1 /3 = 1 /3 = 1<5=1 

(^ (0) " *f )l*&> + # " E l<)<*, ( > (1) l*, ( ?> = (2.67) 

p=i 

Exemple : l'effet Stark 

Si un atome d'hydrogene est soumis a un champ electrique externe, la frequence de la transition 
Lyman-a (Is — > 2p) est deplacee d'une quantite proportionelle au champ electrique, pour des intensites 
de champ inferieures a 10 4 kilovolts/m. 

Ai^cm" 1 ) = -12.8 £ (2.68) 

L'intensite du champ electrique est exprimee en kilovolts/m. Pour des champs plus intenses des effets non- 
lineaires apparaissent. La theorie des perturbations au premier ordre permet d'expliquer ce phenomene 
quantitativement. Le potentiel de perturbation est 

= _ z = r^Y? (2.69) 
et £ est assimile au parametre de perturbation. La correction au premier ordre de l'energie de l'etat s est 

E 1 ^ = (ls\x\ls) = (2.70) 

Le niveau n = 2 est degenere, si Ton choisit d'orienter le champ electrique suivant z, seuls les elements de 
matrice (2s|z|2p 2 ) et (2p^|z|2s) ne sont pas nuls, ils sont egaux a 3 u.a. La diagonalisation de la matrice 
H' 1 ) est immediate et 

E { 2 r) = -3 (2.71) 

done 

Av = {E^P - E { ^)£ = -S.Ou.a. = -12.834 £ (2.72) 
quand est exprime en cm -1 et £ en kvolt/m. 

2.6.2 Perturbation dependant du temps. 

Deux methodes peuvent etre utilisees pour traiter la perturbation dependant du temps: la methode 
de variation des constantes et la methode de variation perturbation. 

Methode de variation des constantes. 

Considerons un systeme quantique se trouvant avant la perturbation it < to) dans l'etat stationnaire 
d'energie auquel correspond: 

|^ 0) (r,*)} = |^(r))exp(-z^/fi) (2.73) 
la fonction |\f, (r, i)} satisfait a l'equation de Schrodinger dependant du temps: 

ih^°\r,t)) = i7(°)(r)|^ 0) (M)} (2-74) 

Supposons maintenant, que la perturbation dependant du temps H^\r,t) soit appliquee a l'instant t = 
et supprimee a l'instant t — r. L'etat du systeme perturbe sera decrit par |\f,(r, i)}. L'indice i n'est pas 



relatit a la situation au systeme apres la perturDation mais, a la situation ongmelle. La tonction a onde 
iVJfj (r, t)} est solution de l'equation de Schrodinger dependant du temps: 

i h^ i (rJ)) = (H(°\r) + H( 1 \rjm i (r,t)) (2.75) 

Quand la perturbation est supprimee le systeme revient a un etat stationnaire et la structure primitive 
en niveaux d'energie est retablie. Pourtant, si avant la perturbation, le systeme occupait le i-ieme niveau 
d'energie, il peut, une fois celle-ci supprimee occuper le /-ieme. Dans ce cas on dit que la perturbation 
H^fojt) induit la transition quantique i — y f. 

Le principe de superposition permet de representer l'etat perturbe |\f, (r, i)} en fonction des etats 
propres de H^: 



*,-(r ,*)) = £ c «WI*i 0) (M)> = ^c ii (i)|$ i (r))exp(-z J B i i/fi) (2.76) 



,(0 

i i 
En substituant 2.76 dans 2.75 on obtient: 



+ £c^)tfM(r,i)|^ 0) (r,*)} (2.77) 



zfi^-c ii (t)|^(r ; t)) + zfi^c ii (t)-|^(r ; t)} = £ c„ <°>(r)|¥< >(r, t)> 
I I I 



I 

En eliminant les termes correspondant a l'equation de Schrodinger non perturbee il reste: 

^Ei.'(i)l* i (0) (M)) = E c «'(^ (1) ( r . i )l*r ) (M)) (2-78) 

En multipliant a gauche par (\Pm^(r, t)\, il vient: 

I 

= ^c ii (i)($( n ° ) (r)|^ (1) (r,i)|$ i (0) (r)}exp( J (i? m - E t )t/h) (2.79) 
I 

Les coefficients Cj/(i) sont alors developpes en serie 

c il (t) = c$\t)+c} l {t) + ... (2.80) 

avec 
done 

|4 0) W = (2.82) 



et 



= 4/ (^ (0) (r)l^ (1) (r,Ol^f ) (r)}exp( J (i? i - Ei)t'/h)dt' (2.83) 



Methode de variation perturbation. 

La methode de variation des constantes ne permet pas de retrouver les equations de la perturbation 
statique, dans la limite ou J?' 1 ) est independant du temps. Dans la methode de variation perturbation, 
l'equation de Schrodinger perturbee s'ecrit: 

[fl^W + AfiWfr.i)] |* i (r,t)> = «fi^|* i (r,t)> (2.84) 



et la ionction propre \Wi\T,t)) est aeveloppee en sene au parametre ae perturoation a: 

CO 

|*i(r, t)) = exp(-^(t)) £ \ k \ X k i (r, t)> (2.85) 

i—k 

ou le facteur de phase (p(t) est lui-meme developpe en serie de A 

CO 

p(*) = ^AV(<) (2-86) 

|Xi°^(r, i)) est choisi tel qu'il soit fonction propre de l'equation de Schrodinger non perturbee: 

|xr ) (r,<)> = |*l° ) (r)> (2.87) 

ce qui implique 

^(o) = E i°) t / h (2. 88 ) 
De plus on impose a |M?i(r, t)) d'etre normee. En substituant 2.85 et 2.86 dans 2.84 il vient: 

(oo \ /oOo\/oo \ 

£A*| X ?(r,*)) eM-M*)) = £ fiA *a^*W £**lx?(M)> exp(-^(t)) 

k=0 / \k=0 J \k=0 / 

+A (x>*^ix?( r >*)>) ex p(-*y(*)) ( 2 - 89 ) 

En conservant les termes de degre homogene en A, on obtient pour le premier ordre: 

(£<°>(r) - E^)\ X l(r,t)) + H^(r,t)\^°\r)) = h^(i)\^°\r)) + ih^\ X }(r,t)) (2.90) 
En multipliant a gauche par ($j c ^(r)|, on trouve : 

h^(t) = (*?\T)\H^T,t)\*l°\T)) - *h(*?\T)\^ X }{T,t)) (2.91) 

En substituant dans 2.90 on obtient: 

(i7(°)(r) - E^)\xl(r,t)) + H^(r,t)\^ \r)) - \^°\r))(^°\r)\H^(r, t)\^%)) = 

- i H\^ \ r ))& \ r )\^xt(rJ)) + 1 h^\x!(r,t)) (2.92) 

Si H 1 est independant du temps, l'equation 2.92 se ramene a celle de la perturbation statique: 

(£<°>(r) " ^ 0) )l*i W> + ^ (1) (r|^ 0) (r)} - |*<°>(r))<*<°>(r)|ff (1) (r)|^ 0) (r)} = (2.93) 

La fonction d'onde perturbee est alors evaluee a l'aide des memes techniques que dans le cas de la 
perturbation independant du temps. 

Exemple: calcul de la polarisabilite dynamique. 

considerons le cas ou le potentiel de perturbation XH^\r,t) = r£(t) est due un champ oscillant de 
frequence u>: 

£{t) = [exp(zwt) + exp(-zwi)] (2.94) 

avec, implicitement 

A = £W (2.95) 



ivietnoae ae variation aes constantes n,n suDstutuant aans 1 equation z.od 

t 

c mi = {<f>° m |r|$°) J [exp(i(E m - Ei + hu)t'/h + exp(i(E m - E> - hu)t'/h)] dt' (2. 



96) 



dans la limite t — > 



et 



c 1 • - ($° |rl$°) 
aH = ^|«|r|$ 



1 



+ 



1 



E m - Ei + hw E m - Ei - hw 
1 1 



0\|2 



+ 



Ei +hw E m ~ Ei - hui 



(2.97) 
(2.98) 



Methode de variation perturbation La fonction d'onde perturbee au premier ordre est ecrite sous 
la forme de la somme de deux composantes: 



|Xi (',<)> = IHWJ^PM) + |*-(r))exp(-«wi) 
En substituant dans 2.92, on obtient pour \$>\) et 

(H° - E° ± Ml*±> - |<D°}(<D°|r|<D°) = 
En developpant sur les fonctions propres de H°: 

|$° \/$° | r |$9\ 



et 



1 ±; 2-> E m - E? ± 

m 

aH = ^|«|r|$°}| ' 



1 



+ 



1 



Ei +hw E m - Ei - hui 



(2.99) 
(2.100) 

(2.101) 
(2.102) 



2.7 Principe variationnel 

La methode variationnelle est une autre technique de calcul de fonctions d'onde approchees dans 
laquelle le principe variationnel est utilisee pour minimiser une fonctionnelle dont la valeur est optimale 
pour la fonction d'onde exacte. 

Les valeurs propres de l'hamiltonien H sont classees par ordre croissant 



Ei+i > Ei 

La plus petite valeur propre, Eq est celle de l'etat fondamental. Si |$) est une fonction d'onde approchee 
qui peut etre choisie normee, alors 

(<f>\H\<f>) > E 

Les fonctions propres de H forment une base complete sur laquelle on developpe |<J>), 



et 



d'ou 



i j i 

(*\H\Q) = E + J2(Ei - E )\(^i\<^)\ 2 > E 



(2.103) 
prochee 

(2.104) 
(2.105) 
(2.106) 
(2.107) 



en enet les lacteurs apparaissant dans la somme sont tous deux positirs. L energie ae 1 etat londamental 
est une borne inferieure a la valeur moyenne de H par rapport a une fonction d'onde approchee. Si |$) 
depend des parametres {(, } les valeurs optimales de ces parametres seront celles pour lesquelles 

Jr($|.ff|$> = (2.108) 

2.7.1 Exemple : 

( 7 1 \ 1 I' 2 

On sait que la fonction d'onde exacte de I'etat fondamental d'un ion hydrogenoi'de est I — J exp (—Zr) , 

7 2 

I'energie de I'etat fondamental est — 4p. Pour retrouver ces resultats a partir du principe variationnel, on 
choisit une fonction d'onde d'essai 



|$)=( — ) exp(-Cr) (2.109) 



7T 



oil £ est le parametre a optimiser 



et 



($|#|$) = t 2 ( 2 ~ Z C (2.110) 



— {<S>\H\9) = C - Z (2.111) 



la valeur optimale de (* est Z , ce qui permet de retrouver la fonction d'onde exacte. 

2.8 Proprietes locales. 

2.8.1 Matrices densite. 

Pour un systeme a N particules de coordonnees d'sepace et de spin xi,X2, . . . ,xjy represente par la 
fonction d'onde \f(xi , x 2 , . . . , xjy) la fonction 

rW(x 1 ,x 2 ,...,x JV ;x / 1 ,x 2 ,...,x^ v ) = *(xi,x 2 ,...,x J v) x **(x'i , A, ■ ■ ■ , *?n) ( 2 - 112 ) 

est appelee matrice densite du systeme totale. Plus generalement, on definira les matrices densite d'ordre 
k par les fonctions: 

rW(xi,x 2 , . . . ; x ft ;x' 1; x' 2; . . .,x£) = jj^-^y J dr k+1 J dr k+2 . . . J drjv^xi , x 2 , . . . , Xjv ) x** (x' 1; x' 2 , . . . , 

(2.113) 

En particulier les matrices densite du second et du premier ordres ont pour expressions: 

r( 2 )(x 1; x 2 ;xi,x 2 ) = N(N-1) jdr 3 J dr 4 . . .J dr N ^( Xl , x 2 , . . . , Xjv ) x **(xi, x' 2; . . . , x^) (2.114) 



et 



r^(xi;x' 1 )=JV J dr 2 J dr 3 ...J dr N ^( Xl , X2 ,..., XN ) x ¥*(xi, *2, • • ( 2 - 115 ) 
La densite electronique est egale aux elements diagonaux de r^^xi ; x'): 

p(r) = r( 1 )(xi;xi) = JV J dr'*(x 1 ,x 2 ,...,x JV ) x * + (xi,x 2 ,...,x^) (2.116) 

ou dr 1 indique que l'integration est effectuee sur les coordonees d'espace et de spin de tous les electrons 
sauf un. 



z.o.z uensite ae propriete. 

Pour la propriete represente par l'operateur A la densite de propriete Pa{^) es t definie par la relation: 

p A (r) = y J dr' <jVi* + (i*)**} (2.117) 

Cas de l'energie cinetique. 

La densite d'energie cinetique A'(r) est definie par la relation: 

K(r) = Y j dr'|^[-Iv 2 ]M/ + (-Iv 2 ]M/r*} (2.118) 

Comme 

V 2 /9(r) = V{(W*)tf + 



alors: 



= (V 2 **)* + **V 2 * + 2V**V* (2.119) 

K(r) = y J dr'WW- ^V 2 p{r) 

= G(r) + £(r) (2.120) 



Chapitre 3 

Symetrie d'espace 



3.1 Groupes de symetrie 

Les noyaux d'ime molecule dans une configuration donnee forment un ensemble de points qui peuvent 
se transformer les uns dans les autres sous Faction d'operations de symetrie. Ainsi, peut-on permuter les 
quatre hydrogenes d'une molecule de methane. 



3.1.1 Groupes ponctuels 

Les operations de symetrie, qui seront notes P,, forment un groupe. Un groupe G est un ensemble fini 
ou infini d'elements pour lesquels on definit une loi de composition interne notee x. Si cette loi satisfait 
les conditions suivante l'ensemble est un groupe. 

1. A chaque couple d'elements P, et Pj correspond un unique element P& 

P k = Pi x Pj (3.1) 

Pk est appele produit de P, et Pj. 

2. la loi de composition interne est associative 

(Pi x Pj) xP k = Pi x (Pj x P„) (3.2) 

3. il existe un element neutre, l'identite E 

Pi x E = E x Pi = Pi (3.3) 

4. chaque element possede un inverse note P _1 

Pi x Pr 1 = P." 1 xP i = E (3.4) 

Si la loi de composition est commutative 

Pi x Pj = Pj x Pi (3.5) 



le groupe est abelien. 

Les elements P, et Pj commutent si 



Pi x Pj = Pj x Pi (3.6) 



L'ordre h du groupe est egal au nombre d'operation du groupe. Soit P, est un elements de G, si toutes 
les puissances de Pj sont des elements distincts de G, l'element P, est d'ordre infini, sinon il existe un 
entier k tel que 

Pf=E (3.7) 
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k est 1 orare ae ae L element neutre est le seul element a orare 1. les orares a un element et ae son 
inverse sont identiques. 

Un sous groupe H est un sous ensemble d'elements de G qui forment un groupe a eux seuls. 
II existe deux types d 'operation de symetrie : 

1. les rotations pures : 

l'identite E 

les rotations C™ de 2mir/n autour d'un axe d'ordre n note C„. L'axe de symetrie 
principal correspond a la plus grande valeur de n. 
Les rotations pures sont des operations propres, c'est a dire qu'elles effectuent des transformations 
congruentes. 

2. les rotations impropres 

l'inversion i 

les plans de symetrie a, ils sont designes par si le plan est perpendiculaire a l'axe 
principal, a v s'il contient cet axe et (Td si c'est le plan bissecteur de Tangle forme par 
deux axes Ci perpendiculaires a l'axe principal, les rotations inversion S n — (rC n qui 
correspondent a la rotation C n suivi d'une reflexion a dans un plan perpendiculaire a 
l'axe, les operations S 1 ™ correspondent a m fois l'operation S n . 
Les operations impropres transforment un objet en son enantiomorphe. 



Symbole de Schoenfiies 


operations de symetrie 




E,C(<f)),(7 v lineaire 


Dcah 


E,C(<f>),C2{a),S(ir + <f>),i, a v ,a h lineaire 


Ci 


E 


a 


E,<7 


Ci 


E,i 


C n 


E, C n 


C n h 




Cnv 


E, C n , u v 


S'2n 


E, C n , S'jn 


D n 


E, C n , C2 


D n h 


E,C n ,C' 2 ,ah 


D n d 


E, C n , C*2, C d 


T 


E, C n , C3 


T h 


E,Cn,C' 3 ,l 


T d 


E, C n , Cg, (Td 





E , C n , C4 


o h 


E, C n , C4, 1 


I 


E, C n ,C 5 


h 


E,C'n,C" 5 ,l 



Tab. 3.1: Signification des symboles de Schoenfiies 

Un groupe d'ordre fini est un groupe fini, un groupe d'ordre infini est dit continu. 

Les differentes operations d'un groupe de symetrie laissent au moins un point invariant, ce point 
est commun aux differents elements de symetrie. Les groupes de symetrie sont appeles groupes ponc- 
tuels pour cette raison. Les groupes ponctuels sont designes par leur symbole de Schoenfiies (Ex : 
Civ, Td, Ih, Dcah,. . .). Les symboles de Schoenfiies sont constitues d'une lettre majuscule C,D,S,0,T,I 
et d'un indice inferieur. Le tableau 3.1 donne la signification des symboles utilises. 

Les groupes ponctuels finis d'ordre n sont isomorphes d'un sous groupe du groupe de permutation S n 
(theoreme de Cay ley). 



o.i.z rjxempie : le groupe U3„ ae ia moiecuie i>n 3 



Les elements de symetrie du groupe permettent de permuter les noyaux des hydrogenes Hi, H2 et H3 
et laissent le noyau de l'azote invariant 1 . La figure 3.1 definit la numerotation des atonies. L'axe C3 est 
l'axe z, les plans a Vl sont definis par l'axe z et les liaisons N-H,. 




Fig. 3.1: Molecule NH 3 

Le groupe Cz v est isomorphe du groupe de permutation S3, aux 3! permutations de 53 correspondent 
les six operations de symetrie. La correspondence entre les operations de symetrie et les permutation des 
atomes d'hydrogene est donne par le tableau 3.2. 



element de symetrie 


operation de symetrie 


permutation 


axe z rotation de 2ir 


E 


(1)(2)(3) 


plan (zNHi) 




(1)(23) 


plan (zNH 2 ) 




(2)(13) 


plan (zNH 3 ) 




(12)(3) 


axe z rotation de 27r/3 


cl 


(123) 


axe z rotation de 47r/3 


c'l 


(132) 



Tab. 3.2: Operations de symetrie et permutation des hydrogenes de NH3 

Le produit de chaque paire d'operations est donne par la table de multiplication du groupe. Le resultat 
de deux transformations congruentes successives est congruent, celui de deux transformations enantio- 
morplie est egalement congruent parceque l'image mirroir de l'image mirroir d'un objet est congruente a 
cet objet. Le type du produit de deux operations de symetri est donne par la regie : 

propre x propre = imropre x imropre = propre 
propre x imropre = imropre x propre = imropre (3-8) 

D'apres la table 3.3, il est clair que le groupe C$ v n'est pas abelien. 

1 une permutation notee [aP^fS) correspond a a — > P — > 7,7 — > S : S — > a 







<J V1 


<JV 2 


<^ 3 


^3 


^3 


E 


E 




< T V 2 


c^ 3 


ci 


eg 






E 


°3 


ci 


&V2 


Cu 3 






ci 


E 


ci 






Cu 3 




ci 


ci 


E 






ci 


ci 






&V2 


E 


eg 


c'i 


ci 


&V 2 


Cu 3 




ci 





Tab. 3.3: Table de multiplication du groupe C3,, 

3.2 Transformation des fonctions d'onde 

Les operations de symetrie transforment les ccordonnees d'espace des particules. Si Ton considere un 
ensemble de N particules reperees par les vecteurs position ri , . . . , rjv, operation P, transforme ces points 
en : 



/ 



PiVl 



(3. 



A chaque operation de symetrie P, on associe un operateur Op t agissant sur l'espace des fonctions des 
variables {r} et defini par l'identite 



Pi f(P i T 1 ,...,P i T N )=f{T 1 ,...,T N ) 



ou / est une fonction quelconque. 
Les operateurss Op sont lineaires 



P (af + hg) = aOpf + bO P g 



et conservent le produit scalaire 



{Opf\Opg) = (f\g) 
ils forment un groupe isomorphe du groupe ponctuel de symetrie, en effet 

Pk = Pi x Pj 

implique 



Les operateurs Op commutent avec H 



done 



Pk = Op, x Pj 
[H,O P ]=0 



P H\^) = HO P \^) = EO P \y) 
si |\?) est fonction propre de H, Op\^) Test egalement pour un niveau non degenere, done 

Op^) = D\V) 

dans cette equation D est une constante multiplicative egale a ±1 si Op|\f) est normee. 



(3.10) 

(3.11) 
(3.12) 

(3.13) 
(3.14) 
(3.15) 
(3.16) 

(3.17) 



o.o xxtijji tistniLciLiuiis 

Pour un niveau /— fois degenere, il existe / fonctions propres lineairement independantes correspon- 
dant a le meme valeur propre E, et dans ce cas 

/ 

Op\^ i ) = Y J D k i{P)\^k) (i = 1,2,...,/) (3.18) 

k=i 

L'equation 3.18 associe une matrice / x / D(P) a chaque operation de symetrie du groupe, en faisant 
agir l'operateur Opi sur 3.18 

/ / / 

0p:0 P \^) = ^£>*i(P)Op.|**) = J2 D ^( p )J2 Dlk ( p, ^ l > 

k=l k=l 1=1 

S S S 

= EE*^'^!*')^^^^)!^) ( 3 - 19 ) 

k=l 1=1 k=l 

done 

D(P' x P) = D(P')D(P) (3-20) 

La matrice correspondant au produit de deux operations de symetrie est egale au produit des matrices 
associees a chacune de ces operations. L'ensemble des matrices D(P,) i = 1, . . . , n) torment une re- 
presentation de dimension / du groupe, et les fonctions une base pour cette representation. Les 
representations sont designees par le symbole T que Ton indice (Ti^T?, . . .) pour distinguer les represen- 
tations les unes des autres. 



3.3.1 Representations reductibles et irreductibles 

Si Ton considere les representations ]?i{Di(P)} et T2{D2(P)}, de dimension fi et on peut 
construire des matrices D(P) de dimension fi + f% 



D(P) = 




(3.21) 



On obtient ainsi une nouvelle representation de dimension fi + f% qui est la somme des representations 
Ti et T?- De meme si toutes les matrices d'une representation T d'une dimension / peuvent etre mises 
sous forme de blocs diagonaux de dimensions fi, fo, . . . par une meme transformation unitaire de la base 
on dit que la representation T est reductible. Une representation irreductible est une representation pour 
laquelle il n'existe aucune transformation de ce type. 



3.3.2 Caractere d'une operation 

La trace de la matrice D(P,) de la representation L est appelee caractere de l'operation P, dans T et 
est notee Xr(Pf)- 

Le caractere d'une operation est invariant par transformation unitaire de la base, en effet si U est une 
matrice unitaire 

trUtDU = trD (3.22) 

si 

r = L 1 + L 2 alors Xr (P) = Xr, (P) + Xr 2 (P) (3.23) 

Les representations irreductibles sont designees par leur symbole de Mulliken compose d'une lettre ma- 
juscule et d'indices superieurs et inferieurs. La correspondance entre les caracteres des operations et les 
symboles de Mulliken est indiquee dans le tableau 3.4 Pour chaque groupe les caracteres des operations 



sont aonnes pour criaque representation irreductiDle par la taDle ae caractere au groupe. La premiere 
ligne de la table de caracteres indique les operations du groupe, la premiere colonne les representations 
irreductibles, les deux dernieres colonnes indiquent la symetrie des fonction de x, y et z des premiers et 
second degres centrees au point invariant. La representation irreductible dont tous les caracteres sont 
egaux a 1 est appelee representation totalement symetrique. La table de caracteres du groupe C^ v est 
reproduite dans le tableau 3.5. 



symbole de Mulliken 


carcteres des operations 




771 
t, 




C2{-LC n ) C 2 (±C n ) C 2 (1C„) 




lettre 










A 


1 


1 






B 




"I 






Hi 


I 








77 

r 


3 








indice inferieur 










1 






1 




2 






1 




3 
9 






i 
i 


1 


u 








-1 


indice superieur 
/ 








1 


// 








-1 


+ 








1 
-1 



Tab. 3.4: Correspondance des symboles de Mulliken et des caracteres pour les groupes finis 
La representation totalement symetrique du groupe C$ v est la representation A\. 



C*3u 


E 


2C 3 


3d"„ 








1 


1 


1 


z 


x' 2 + y' 2 + z' 2 


M 


1 


1 


-1 






E 


2 


-1 





( X : v) 


(x 2 - y 2 ,xy), (xz,yz) 



Tab. 3.5: Table de caractere du groupe Cz v 



3.3.3 Reduction des representations 

Une representation reductible L est la somme de representations irreductibles L, 

L = ^m i L i (3.24) 

i 

le caractere de l'operation P dans la representationr est done 

X (P) = "£miXr,(P) (3-25) 

i 

inversement les coefficients mi peuvent etre calcules en appliquant la relation 

m i = \Y.*h(P)x{P) (3-26) 
p 



li est 1 ordre cm groupe et la somme est enectuee sur toutes les operations au groupe. 

Le produit direct de representation est une application importante de la reduction de representations. 
Si les fonctions {<f>x, <f>2, ■ ■ ■ , 4>j} forment une base de la representation Li de dimension / et les fonctions 
{fiifii ■ ■ ■ , fg} une base de la representation L 2 de dimension g, les fonctions 4>kfm {k — 1, f;m — 
1, </) forment une base de dimension fg de la representation produit direct 

ri «>r 2 

Les caracteres de la representation Li ® L 2 sont egaux aux produits des caracteres des representations 
d'origine 

Xr 1 «r 2 (P) = Xr 1 (P)xr 2 (P) (3.27) 

Ce resultat est utilise pour construire, en appliquant l'equation 3.26 des tables qui donnent les reductions 
de representation des produit directs des representations irreductibles d'un groupe ponctuel. La reduction 
de representation des produits directs des representations irreductibles du groupe C$ v est donnee dans le 
tableau 3.6. 





Ax 


A-2 


E 




A x 


A 2 


E 


A 2 


A 2 


Ax 


E 


E 


E 


E 


Ax + A 2 + E 



Tab. 3.6: Reduction des produits directs des representations du groupe C$ v 



3.4 Utilisation de la symetrie moleculaire 

L'operateur hamiltonien d'un atome ou d'une molecule est invariant par permutation des noyaux done 
par rapport aux operations du groupe ponctuel, ses fonctions propres forment des bases pour toutes les 
representations irreductibles du groupe. Dans le cas d'etats degeneres, la degenerescence est egale a la 
dimension de la representation correspondante. 

3.4.1 Evaluation d'integrales 

Une propriete importante concerne la valeur d'une integrale de la forme 

</|G|<7> 

Si f,g et G sont des vecteurs de base des representations Li,L2 et L3, l'integrale (f\G\g) est non nulle si 
le produit direct 

Li ® L 2 ® L 3 

contient la representation totalement symetrique du groupe. 

Cette propriete permet de deduire les regies de selection en spectroscopic, en effet, la probabilite de 
transition entre les etats i et j est proportionelle a l'integrale 

dans ce cas G est l'operateur qui couple 1 'hamiltonien non-perturbe du systeme au champ externe qui in- 
duit la transition. Dans le cas des transitions photoniques cet operateur est l'operateur moment dipolaire. 

et sont fonctions propres de H et done fonctions de base de representations irreductibles L, et 
Tj du groupe. L'operateur G peut toujours etre associe a une representation reductible ou irreductible 
du groupe, Tq si le produit direct 

Li ® L G ® Tj 



ne contient pas la representation totalement symtrique au groupe 



(*i|G|*j> = (3.28) 

et la transition est dite interdite. 

3.4.2 Construction de fonctions de symetrie 

Pour construire des fonctions de symetrie a partir de fonctions de base, il faut d'abord construire la 
base, puis proceder a la reduction de la representation decrite par la base et enfin trouver la transformation 
unitaire correspondante. Pour construire la base on choisit une ou plusieurs fonctions sur lesquelles on 
fait agir les operations du groupe. Dans le cas de la molecule NH3, par exemple, a partir d'une fonction 
ls(l) centree sur Hi les differentes operations engendrent les fonctions suivantes : 



B 1,(1) 


= Ml) 


O c ,l,(l) 


= 1,(2) 


O c As(l) 


= 1,(3) 




= Mi) 




= 1,(2) 




= 1,(3) 




(3.29) 

La base est de dimension trois et les matrices de transformation sont : 

/ 1 

D(E) = 1 I D{Cl) 
\ 1 

1 

D{Cl) = I 1 I D(a Vl ) 
1 

1 

D{a V2 ) = I 1 ) D{a Vs ) = \ 1 | (3.30) 

1 

Les caracteres des differentes operations sont 

X (E) = 3, X(C 3 ) = 0, X(°v) = l (3.31) 
En appliquant l'equation 3.26, on trouve : 

m Al = i(lx3 + 2xlx0 + 3xlxl) = l 
m A2 = i(lx3 + 2xlx0 + 3x-lxl)=0 

m E = -(2 x3 + 2x-lx0 + 3x0xl) = l (3.32) 
6 

done 

T = A 1 +E (3.33) 
Pour trouver les fonctions de symetrie (f)^ , on utilise la relation 

4 n =E4V)Oi^ (3-34) 



ou ip est une tonction de base generatnce, et les L)\^ \f) sont les elements de matrices dermis par 1 equation 
3.18 pour la representation irreductible T Pour le groupe C$ v ces matrices sont donnees dans le tableau 
3.7. On obtient par cette methode 



<^ lJ = 2(1«(1) + 1«(2) + 1«(3)) 

<f>ii 2) = o 

<j>[f = 21«(1)-1«(2)-1«(3) 

4,W = v / 3(ls(2) - 1«(3)) (3.35) 



A 1 1 1 1 1 1 1 

A 2 1 1 1-1-1 -1 




Tab. 3.7: Representations irreductibles de C$ v 



Chapitre 4 



Separation des mouvements 
electroniques et nucleaires 



4.1 Position du probleme 

Une molecule est constitute de deux types de particules : 

- les noyaux de masse de charge Z^e, qui peuvent etre des fermions ou des bosons 

- les electrons de masse m e , de charge — e qui sont des fermions. 

En unites atomiques, l'hamiltonien correspond ant a un tel systeme a pour expression : 

? 2 i — — 1 - — 7- — Z A Z B 



h= -\Y.jn - sEv? + E^-EEfr+E 



A A i i<j 13 i A lA A<B Aa (4.1) 



Tel V{B.,T) 



L'operateur H classique depend uniquement des variables d'espace : les coordonnees {Rji} des noyaux 
et {r;} des electrons. Le caractere fermionique ou bosonique des particules n'apparait pas dans l'operateur. 

Le premier probleme a resoudre va etre d'essayer de separer le mouvement des noyaux de celui des elec- 
trons. Ceci est d'autant plus justifie que ces particules ont des masses tres differentes et qu'elles peuvent 
avoir des spins differents. H peut etre decompose en une contribution cinetique purement nucleaire Tjv 
et un hamiltonien electronique : 

H el =f el + V(R,r) (4.2) 

qui est un operateur hermitique parametrique en R. L'ensemble de ses fonctions propres {i/>(R ; r,E)} - E 
designe l'ensemble des variables de spin electronique a, - forme une base complete sur laquelle on pourra 
developper la fonction d'onde totale du systeme |\?) 

|¥) = £*„.(R)|lMR,r,£)) (4.3) 

n 1 

ou les $(R) sont des coefficients parametriques. En substituant dans l'equation de Schrodinger totale on 
obtient : 

J2{T N + Hel}$n' (R)|VVi'(R, r,E)> = E ( R ) 1^"' ( R , r , S )) ( 4 - 4 ) 
n' n 1 

comme : 

H el \i> n , (R, r,E)> = E n , \i> n , (R, r,E)> (4.5) 

^{f J v + J B n ,}* n ,(R)|V n .(R,r,E)> = J B^* n .(R)|V n -(R,r,E)> (4.6) 
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en multipliant a gaucne par [ip n (ti.,T,2j)\ on oDtient : 

{f N + U n (R) - E} \i> n (R, r,E)> = - £ C nn /(R, P)$ n , (R) (4.7) 

avec : 

U n (R) = E n (K) + A n (R) (4.8) 

A n (R) = (^ n |T JV |^n) (4.9) 

et 

C nn ,(R,P) = (i> n \T N \i> n ,) + Y / -±-((i> n \P A \i> n ,) (4.10) 

A MA 

A cause du second membre on doit connaitre le spectre complet des etats propres de H e i . 

4.1.1 Approximation adiabatique 

Dans l'approximation adiabatique on suppose que 

£ C nn .(R,P)*„.(R) = (4.11) 

On a done a resoudre l'equation de Schrodinger electronique pour differentes configurations nucleaires, 
puis l'equation nucleaire : 

{f N + U n {B.)}\*) = E\*) (4.12) 
La fonction d'onde totale s'ecrit sous la forme d'un produit : 

i*> = \®m (4.i3) 

|$) peut maintenant etre identifiee a la fonction d'onde nucleaire et \ip) a la fonction d'onde electronique. 
Cette approximation est excellente si : 

m separation des niveaux de vibration v 

X : : : f» 10" 7 (4-14) 

M separation des niveaux electroniques 

4.1.2 Approximation de Born-Oppenheimer 

Dans l'approximation de Born Oppenheimer on neglige A„(R) et l'equation de Schrodinger nucleaire 
devient : 

{f N + E n (R)}\$) = E\$) (4.15) 
La surface definie par £" n (R) est appelee surface de Born-Oppenheimer. 



Chapitre 5 

Une molecule monoelectronique H2 



L'ion moleculaire est la molecule la plus simple. Dans l'approximation adiabatique, il existe des 
solutions analytiques exactes a l'equation de Schrodinger electronique. De plus, dans les cas limites 
corespondant aux distances internucleaires nulles ou infinies, ces solutions s'expriment simplement a 
partir des orbitales atomiques des ions hydrogenoi'des. La nomenclature utilisee pour nommer les orbitales 
moleculaires de a d'abord ete etendue aux molecules diatomiques, puis a toutes les molecules. 



5.1 Fonction d'onde exacte 

L'hamiltonien electronique de l'ion moeculaire est 



H = 



1 1 



n R 



(5.1) 



ou R designe la distance internucleaire, r a et rj, les distances entre l'electron et les noyaux A et B. pour 
resoudre l'equation de Schrodinger, on se place en coordonnees elliptiques ^,r},f. 



T] = 



r-g + n 

R ' 
r a - r b 

R ' 



1 < £ < co 

-1 < v < 1 



(5.2) 



if est Tangle azimuthal autour de l'axe internucleaire, son domaine de variation est 

0<<p<2ir 

Dans le repere cartesien, les noyaux A et B ont pour coordonnees respectives (0, 0, —R/2) et (0, 0, -R/2). 
En coordonnees elliptiques l'hamiltonien s'ecrit 1 



H = 



d 



(e - 1) 



d_ 
dt] 



dr) 



e-v 2 



d 2 



+ 2ZR-§(e-V 2 ) 
(5.3) ' 



En coordonnees elliptiques, la fonction d'onde electronique s'exprime sous la forme d'un produit 

Vtf,iJ,p) = G(0i7fo)*te) (5.4) 
Les fonctions G(£), H(rj) et $[y>) sont solutions des equations differentielles 



d 

dk L 



(e - 1) 



d 

dl\ 



+ 



X 2 C 2 + 2R£ 



G(0 = 



(5.5) 



^^L'expression du laplacien dans differents systemes de coordonnees est donne dans l'annexe A. 
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drj 



\ 2 2 
AT) 



(1 - ^) 



fffo) = 



— $(<p) +m 2 $(v?) = 
ay- 4 



dans cette equation ji est la constante de separation, et 



(5.6) 

(5.7) 
(5.8) 



(5. 



La composante suivant l'axe internucleaire de l'operateur moment angulaire est L z = — i-g^, comme L z 
et H commutent, $(y>) doit etre fonction propre de L, 

L z 4>(<p)=m<l>(<p), (m = 0,±l,±2,...) (5.10) 

et 



$(<£>) = exp (imip) 



(5.11) 



Par analogie avec la nomenclature des orbitales atomiques (s,p, d, f, . . .) les orbitales moleculaires de 
seront designees par les lettres grecques <r, ir, 5, <f>, . . . pour m = 0, ±1, ±2, ±3, . . .. Les formes analytiques 
des fonctions G(£) et -ff(?7) sont beaucoup plus compliquees 



G(0 = exp (-rf)tf + lf(£ 2 - l)l m l/ 2 6(0 



(7 = m — 1 

P 



et 



n = 

Les coefficients C n satisfont la relation de recurrence 

C n — l = U n C n ~\~ v n C n 4. 1 



(5.12) 



(5.13) 



(5.14) 



u n +i = 

v„+i = 

1 = 

V = 



Hirf) a pour expression 



[2n(n + 1) + 2(2p - a + l)(n + 1) - q] 
[n(n - 1) - n(2cr + |m| - 1) + a(cr + |m|)] 
(n + 2)(n\m\ + 2) 

[n(n - 1) - n(2(7 + |m| - 1) + a(a + |m|)] 
-p 2 + 2pcr + cr(\m\ + 1) + v 
\m\{\m\ + 



H(r,)= £ a,pH(,,) 

;=|m| 



(5.15) 



(5.16) 



ou les pj m][ sont les fonctions associees de Legendre d'ordre / et de degre \m\. La relation de recurrence 
des coefficients a; est 



(5.17) 



"1 + 2 



(2/ + l)(2/ + 3) 



{I - \m\ + - \m\ + 2) 

'(1 + 3- |m|)(/ + 3+|m|) (1 + 2 



i|)(i + 2+|m|) /i + (/ + 2)(/ + 3) 



1+2 - 



(2/ + 5)(2/ + 7) (2/ + 3)(2/ + 5) 

(21 + l)(2i + 3)(i + |m| + 3)(i + |m| + 4) 



(5.18) 



(I - \m\ + l)(i - |m| + 2)(2i + 7)(2i + 9) 

D'apres l'equation 5.17, il existe deux series possibles de coefficients a;, 1'une correspondant a / pair, 
l'autre a / impair. 

Pour R — 0, les fonctions d'onde de l'ion moleculaire et de l'ion He + sont identiques. Les fonctions 
de l'atome uni seront utilisees pour nommer les differents etats de la molecule. 

5.1.1 Symetrie des orbitales moleculaires 

Le groupe de symetrie de l'ion moleculaire Hjj" est le groupe ponctuel D^^. Les operations du groupe 
ainsi que leur action sur les coordonnees elliptiques sont presentees dans le tableau 5.1 Les representations 



operation 


description 




Rt) 


Rip 


E 


identite 


n 




f 


C(±<t>) 


rotation de +<j> autour de l'axe 




V 


<p±(f> 




internucleaire 








C(tt) 


rotation de 7r autour de l'axe 






if + TT 




internucleaire 








CM 


rotation de 7r autour d'un axe 




-n 


a — if 




du plan xy et d'azimuth a/2 








% 


inversion 




-n 


if + TT 


S(tT ± <f>) 


S(tt + 4>) = % C{4>), S(n - <f>) = iC(-<f>) 


n 


-v 


If + TT ± <f> 


<*h 


plan mirroir xy, cr^ — iC(tt) 


n 


-v 


f 


a v 


plan mirroir contenant l'axe C(+<j>) 


n 


V 


a + tt — if 



Tab. 5.1: Operations de symetrie du groupe ponctuel 

irreductibles et les caracteres des operations du groupes sont presentes dans le tableau 5.2 La table 5.1 
indique que les differentes operations du groupe laissent £ invariant, mais par contre que rj change de 
signe sous l'effet de l'inversion, des rotations C' a (Tr), des rotations inversion S(TT — +<j)) et des plan mirroirs 
cr v . La symetrie de la fonction d'onde est determinee par les fonctions Hirf) et <&(if). Dans le cas ou les 
coefficients du developpement de Hirf) non nuls correspondent aux valeurs paires de / 



H{~n)= y: a^ H (-»»)= i:(-iy +|m| p i |mi c?)=(-i) |mi ^('?) 

/=|m| /=|m| 

dans le cas ou ces coefficients sont impairs 

OO CO 

H{-t,) = £ a l P} m \-r ] )= ^(-ly+MpM^^M+i^) 

/=|m| /=|m| 



(5.19) 



(5.20) 
(5.21) 



D'autre part, 

$(ip + tt) = exp (im<p) exp (mnr) = (-l) m $(tp) 

Les fonctions g, c'est a dire symetriques par rapport a l'inversion sont sont celles qui correspondent 
a une meme parite de m et des coefficients du developpement de Hirf). Les fonctions u antisymetriques 
correspondent a des parites differentes 2 . Les orbitales moleculaires de sont A ou E suivant que m est 

2 Les indices g et u sont les abreviations de gerade et ungerade. 



oil non mil, g ou u suivant que la tonction est symetrique ou antisymetrique par rapport a 1 inversion. 
Les symboles a, ir, . . . qui caracterisent la valeur de |m| sont utilises a la place des symboles de Mulliken 
A et E. 



representation 


E 


2C(tf) 


C£(tt) 


I 


S(tt ± 0) 


cr„ 


A lg £+ 


1 


1 


1 


1 


1 


1 


Am £+ 


1 


1 


-1 


-1 


-1 


1 


A 2g S" 


1 


1 


-1 


1 


-1 


-1 


Mu s- 


1 


1 


1 


-1 


-1 


-1 


£-i g n g 


2 


2 cos <f> 





1 


2 cos <f> 





-^lu n u 


2 


2 COS (2) 





-1 


-2 cos ^ 







2 


2 cos 2</) 





1 


2 cos 2^ 





a 


2 




o 


_1 


2 cos 2<h 


o 




2 


2 cos 3</) 





1 


2 cos 3</) 





E 3u 


Z 


2 cos 3<p 





-1 


-2 cos 3<p 







2 


2 cos \<j> 





1 


2 cos \<j> 







2 


2 cos 





-1 


-2 cos 






Tab. 5.2: Representations irreductibles et caracteres des operations du groupe ponctuel 

A chaque orbitale moleculaire correspond un etat de l'atome uni qui est utilise pour designer l'etat 
de . Ainsi a-t-on pour les dix premiers etats 

lscr g: 2scr g: 3scr g: 2pcr U: 3pcr U: 4pcr u ,3dpcr g , 4fpcr u ,3pTT U: 3dTT g (5.22) 



atome uni 


atonies separes 


lSCTg 


a(U) 


2sa g 


a(2s) 


Zscr g 


cr(3s) 


2pa u 


a*(ls) 


3pa u 


(7* (2s) 




(7* (3s) 


3d(7 g 


a(2p) 




a*(2p) 


2pir u 


7v(2p) 


3dTT g 


n*(2p) 



Tab. 5.3: Correlations entre les etats de l'atome uni et les etats des atomes separes de l'ion Hjj" 



5.1.2 Orbitales liantes et antiliantes 

Le comportement des orbitales molecul aires pour R^oo fournit des renseignements sur leur caractere 
liant ou antiliant. Par exemple les etats ls<T g et 2pa u correspondent tous les deux a l'etat Is des atomes 
separes comme le montrent les courbes de l'energie electronique en fonction de la distance internucleaire 
(figure 5.1). 

En tenant compte de la symetrie par rapport a 1 'inversion on peut ecrire 

lscr g Si ls a + lsj, si cr 9 (ls) = c(ls) 

2pa u si ls a — lsj, si o"„(ls) = cr*(ls) (5.23) 



Les symDoles cr g [is), cr[Ls), cr u {Ls) et cr [is) sont ceux utilises pour designer les orDitales moleculaires 
en fonction des etats des atonies separes. L'asterisque indique que l'orbitale est antiliante. Une orbitale 
antiliante est une orbitale pour laquelle le plan bissecteur cr^ est une surface nodale, c'est a dire une 
surface ou la fonction d'onde s'annule pour changer de signe. Le plan cr^ est le plan r\ = 0. Si Hirf) est 
une fonction impaire cr^ est un plan nodal, l'orbitale est antiliante, si Hirf) est paire cr^ n'est pas un plan 
nodal et l'orbitale est liante. 

Pour les orbitales liantes la derivee de l'energie electronique par rapport a la distance internucleaire 
est positive au voisinage de R = 0. 



4fo-„ 4po-„ 




0.0 2.0 4.0 6.0 8.0 



T» t.. „ \ 

Fig. 5.1: Energie des etats purement electroniques de H^" 



5.2 Fonctions d'onde variationnelles 

II n'existe pas de fonctions d'onde analytiques exactes pour des systemes autres que les ions molecu- 
laires monoelectroniques. II est done necessaire dans le cas de systemes polyelectroniques d'utiliser des 
fonctions d'onde approchees. Les performances de diverses techniques d'approximation peuvent etre tes- 
tees sur l'ion moleculaire H^" dont la fonction d'onde exacte constitue une reference. 

5.2.1 La fonction de Guillemin et Zener 

La fonction d'onde approchee la plus simple pour l'etat fondamental de Hjj" qui tend vers la fonction 
d'onde exacte de l'atome uni (R = 0) et des atomes separes, est la fonction de Guillemin et Zener a un 
parametre variationnel. L'addition d'un second parametre variationnel ameliore de facon tres significative 



la description de 1 etat tondamental aux distances mtermediaires. oette tonction a pour expression en 
coordonnees elliptiques : 

^ g = N exp (-«{) cosh (3rj (5-24) 

oil a et /3 sont les parametres variationnels. Le facteur de normalisation N est calcule en utilisant les 
integral es A n et B n fournies dans l'annexe B. 



»r-2 

N = 



A 2 (2a) [B (2p) + B o (-20) + 4] - A (2a) 



B 2 (2p) + B 2 (-2p) + j 



(5.25) 



16 

Les parametres variationnels sont fonction de la distance internucleaire R, la valeur optimale est 

( R=2.0 
1 a =1.36 
( /3 = 0.45 

L'energie electronique calculee avec ces parametres est superieure de 2. 10 -4 u.a. a l'energie exacte. 

5.2.2 La fonction de James 

La fonction de James est une fonction a deux parametres variationnels 

* g = iVexp(-0(l+^ 2 ) (5.26) 

qui donne aux distances intermediaries et pour l'atome uni des resultats de qualite comparable a la 
fonction de Guillemin et Zener. La fonction de James n'a pas la forme asymptotique correcte pour R 
grand. 

On remarquera que la dependance en r\ de la fonction de James correspond aux deux premiers termes 
du developpement en serie du cosinus hyperbolique. Les parametres optimaux de la fonction de James 
sont : 

( R = 2.0 
< a =1.30 
( /? = 0.44475 

5.2.3 Les fonctions LCAO 

La fonction de Guillemin et Zener a un seul parametre a — /3 correspond a la combinaison de deux 
orbitales Is, d'exposant (* = 2a/ R : 

V g = N (exp (-C(r - R^) + exp (~C(r - R B )) (5.27) 

ou r, Tlji et designent les vecteurs position de l'electron et des noyaux. La valeur optimale de l'exposant 
C = 1.228 permet de trouver une energie 1.5 10 -2 u.a. superieure a l'energie exacte et de prevoir la distance 
internucleaire correcte. 

Si Ton supprime le parametre variationnel, en fixant £ — 1.0, la distance d'equilibre calculee augmente 
de 0.5 u.a. et l'energie variationnelle de 2.22 10 -2 u.a. 

5.2.4 Les developpements monocentriques 

Une autre methode de construction de fonction d'onde approchee est le developpement sur une base 
finie et done incomplete de fonctions centrees au milieu de la liaison. Avec des fonctions de base de Slater 

Sm m = (2Cr +1 ' 2 [(2iV)!]" 1/2 r-^nfi) (5.28) 

L'etat fondamental est de symetrie <r g , le developpement ne contient que des fonctions de base pour 
lesquelles / est pair et m — 0. On trouve les resultats reportes dans le tableau 5.4 qui montrent que la 



convergence est assez rapide. Le comportement asymptotique ae la convergence en ionction au plus grana 
nombre quantique orbitale utilise dans le developpement est de la forme : 



(5.29) 



developpement 


energie electronique 


convegence e; 


4s 


-1.01850 




4s4d 


-1.08367 


0.06517 


AsAd2g 


-1.09640 


0.01273 


4s4d2g2i 


-1.09992 


0.00354 


4s4d2g2i2k 


-1.10121 


0.00127 


exact 


-1.10263 





Tab. 5.4: Energie electronique de calculee par des developpements monocentriques. 



5.3 Energie locale et potentiel de fluctuation 




Fig. 5.2: Potentiel de fluctuation de la fonction LCAO (C = 1) de H^" 



On peut rendre compte de la qualite des fonctions d'onde approchees en etudiant les fluctuations de 
l'energie locale. Pour une fonction d'onde approchee, l'energie locale est une fonction des coordonnees 
instantanees des particules 



E(r) 



*(r) 



(5.30) 



fom une ionction d onde exacte, 1 equation de Scnrodiger implique que 1 energie locale est constante. Le 
potentiel de fluctuation est la difference entre l'hamiltonien exact et l'energie locale. 

Vfluet {T) = H--^±± (5.31) 

La figure 5.2 montre que la fonction LCAO simple est tres mauvaise pour les faibles valeurs de £ et en 
particulier au voisinage des noyaux ou elle diverge. 



Chapitre 6 



Le spin electronique dans les 
molecules polyelectroniques 

6.1 Theorie du spin de Pauli 

6.1.1 Les fonctions de spin 

Dans la theorie classique du spin electronique on introduit une variable de spin a qui ne peut prendre 
que les valeurs discretes 

* = ±\ (6.1) 
et deux fonctions de spin a(cr) et /3(c), qui forment une base complete et ont pour valeur : 

a(i) = l a(-i) = 

A 

Ces fonctions sont orthonormees : 



/?(„) = o /?(-!) = i 



/ a(a)a{cr)da = a 2 {cr) = 2 + l 2 = 1 (6.2) 
I (3{a)p{a)da = £ /?» = l 2 + 2 = 1 (6.3) 
I a(cr)P(cr)dcr = a(a)(3(a) = 0x1 + 1x0 = (6.4) 

et sont fonctions propres des operateurs s 2 et s z 

3 1 
s 2 a(cr) = — Ti a(cr) s z a(cr) = — Tia(cr) 

6.1.2 Representation matricielle de l'operateur moment cinetique de spin 

L'operateur moment cinetique de spin note s 



S 



(6.5) 
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possede trois composantes qui satistont les relations de commutation : 

[s x , s y] — l h g z (6-6) 

[s y , s z ] = ihs x (6.7) 

s x] = ihsy (6.8) 

on introduit egalement les operateurs de creation et d' annihilation 



s + — s x + is y s+a(cr) = s_|_/?(o") = Tia(a) 

s_ = s x — iSy s_a(cr) = h/3(a) S-/3(a) = 



la representation matricielle de ces operateurs est 



, _ 3^ / 1 A _ ft / 1 

~ 4 I 1 J Sz ~ 2 V -1 

o i A , / o o 



o o J \ 1 o 

1 \ ft / -i 



Sx - 2 V 1 I " y - 2 V » 



(6. 



6.2 Principe d'antisymetrie et ses consequences. 

6.2.1 Principe d'antisymetrie. 

Pour un systeme a deux fermions la fonction d'onde totale (fonction des variables d'espace et de spin) 
est antisymetrique par rapport a la permutation des deux particules : 

*(ri,r2,eri,er 2 ) = -*(r 2; r i; a 2 , vi) (6.10) 

H et s n'agissent pas sur les memes coordonnees done 

[ J ff,s] = (6.11) 

La generalisation de ce postulat aux systemes polyelectroniques de dimension superieure constitue le 
principe d'antisymetrie ou principe de Pauli : 

La fonction d'onde d'un systeme a n particules est antisymetrique par rapport a la permutation de 
deux particules quelconques. Si Ton designe par fi,- l'ensemble des coordonnees d'espace et de spin de la 
particule i, alors 

¥(fli, Q 2 , ■ ■ ■ , Oi, • • • , fy, • • • , ^) = "2, ....fi.; fi„) (6.12) 

Le principe de Pauli est un principe d'exclusion, en effet, si deux particules ont les memes coordonnees 
d'espace et de spin fi,- = fij, alors 

*(I2l, fi 2; . . . , fij, . . . , fij, . . . , fi„) = -*(J2l, S2 2 , ■■■■fi; fi/ ^) 

= ¥(fi 1 ,fi 2 ,...,fi i ,...,fi i ,...,fi„) (6.13) 

ce qui ne peut etre satisfait que si 

*(J2i,J2 2 ,...,S2i,...,S2i,...,S2 n ) = (6.14) 

Deux fermions de meme spin ne peuvent pas occuper la meme position. Le spin d'une particule cree 
autour d'elle une zone d'exclusion pour les particules de meme spini appele trou de Fermi. 



o.z.z wperaxeurs ae spin total 

Pour im systeme a n particules on definit le moment de spin total en sommant sur toutes les particules 
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(6.15) 



on definit egalement les operateurs S+ et S- 

Ces differents operateurs commutent avec H et les fonctions propres de H doivent l'etre egalement 
de S 2 ainsi que de l'une des composantes de S, on choisit en general la composante suivant z. 

A partir des operateurs S x et S y on construit les operateurs 



S- 



(6.16) 



En designant par 6(5', Ms) les fonctions propres de S 2 , Taction des different operateurs s'ecrit 



s 2 e(s, m s ) = s(s + i)n 2 e(s, m s ) 



s z e(s, m s ) = M s ne(s, m s ) (6.17) 



s+e(s,M s ) = y/s(s + 1) - M S (M S + i)ne(s, m s + 1) 

S.Q(S,M S ) = y/S(S + 1) - Ms(M s - l)fi0(5, M s - 1) (6.18) 

L'operateur S 2 ne peut pas etre ecrit sous la forme d'une somme d'operateurs monoelectroniques, mais 
ilpeut etre exprime a partir des operateurs lineaires 5^ , S- et S+ . 

S 2 = S-S++S 2 + TiS z (6.19) 
S 2 = S+S- +S 2 -TiS 2 (6.20) 



6.2.3 Cas des systemes a deux electrons 

Dans les systemes a deux electrons ont peut former a partir des quatre fonctions de spin possibles: 

a(«7i)a(«7 2 ); a(«7i)/?(<7 2 ); p{^)a{^); /%i)/% 2 ) (6.21) 
une fonction normee antisymetrique 



6(0, 0) = ^=[a[a{)p[a 2 ) - /%i)a(<r 2 )) (6.22) 



et trois fonctions symetriques 



0(1,1) = a(ai)a(a 2 ) 

6(1,0) = ^a[v{)p[v 2 ) + p[v{]a[v 2 )) 

6(1,-1) = /%i)/% 2 ) (6.23) 

Les fonctions d'espace <£>(ri,r 2 ) solutions de l'equation de Schrodinger peuvent etre symetriques ou an- 
tisymetriques par rapport a la permutation des coordonnees electroniques ri et r 2 . Les fonctions d'onde 
^(ri, r 2 , <7i, cr 2 ) satisfaisant le principe d'antisymetrisation sont construites en multipliant une fonction 
d'espace par une fonction de spin de symetrie opposee. 



^as aes sysiemes a pius ae aeux electrons 

Pour un systeme a trois electrons il peut exister des fonctions propres de H symetriques par rapport 
a la permutation des coordonnees d'espace : 



<p(ri,T 2 ,T 3 ) = <p(ti,t 3 ,t 2 ) = v?(r 3 ,r 2 ,ri) = y>(r 2 , ri, r 3 ) 
Pour obtenir une fonction d'onde totale 

*(ri,r2,r 3 ,cri,cr 2 ,cr 3 ) = v(ri,r 2 ,r 3 )6((7i,(72,o- 3 ) 
qui satisfasse au principe de Pauli, la fonction devra avoir la propriete suivante : 

6(0-1,0-3,0-2) — e(o- 3 ,o- 2 ,(7i) = e(0- 2 ,0-i,o- 3 ) = -0(0-1, 0- 2 , 0-3) 



(6.24) 
(6.25) 
(6.26) 



II est impossible de construire une telle fonction a partir des fonctions de spin a(cr) et fiia), aussi la 
fonction d'espace totalement symetrique n'a aucune signification physique. 

Cet exemple montre qu'une autre methode doit etre trouvee si Ton veut pouvoir traiter des systemes 
polyelectroniques . 

6.2.5 Spinorbitales et determinants de Slater 

La methode des determinants de Slater permet de resoudre ce probleme dans le cas le plus general. 

Spinorbitale 

Une spinorbitale est une fonction monoelectronique des variables d'espace et de spin. De maniere tres 
generale cette fonction a pour expression : 



<f>(x) = ip + {r)a(a) + ip_(r)/3(a) 



(6.27) 



Les parties d'espace <p+(r) et <fi-(r) sont appelees orbitales. Dans la plupart des applications on utilise 
une version simplifiee de la formule precedente : 



4>(x) 



<p(r)a(a) 

ou 
<p{r)fHv) 



(6.28) 



ou une meme orbitale <p(r) permet de construire deux spinorbitales orthogonales. On trouve souvent la 
notation : 

<t>{x) = <fi(r)a(a) <j>{x) = f{v)P{cr) (6.29) 
qui permet de distinguer les deux spinorbitales. si <j> et <j> sont normees alors 



J <f>*(x)<f>(x)dx = j <f>*(x)<f>(x)dx = 



(6.30) 



Determinant de Slater. 



On appelle determinant de Slater d'ordre TV un determinant d'ordre N formee sur N spinorbitales 
distinctes. 

^(l) <p k2 (i) ... 4> kN (i) 

4> kl {2) <f> k2 (2)l ... 4> kN {2) 



D 



K 



1 



<t> kl {N) <j> k2 {N) 



(6.31) 



le lacteur est un lacteur ae normalisation valaDle si les spmorDitales sont elles-meme normees. A la 
permutation des coordonnees d'espace et de spin de deux electrons correspond la permutation des deux 
lignes correspondantes du determinant. Les determinants de Slater satisfont le principe d'antisymetrie, en 
effet un determinant change de signe lorsque Ton permute deux lignes ou deux colonnes. Cette propriete 
est egalement valable pour une combinaison lineaire de determinants. 
Differentes notations sont utilisees pour les determinants de Slater : 

Dk = -L=-{*i(l),*2(2),---,Mtf)} 
Viv ! 

= |Ai(l), A 2 (2), • ..,k N (N)\ 

= Ak 1 (l)k 2 (2)...k N (N) (6.32) 
Cette derniere expression utilise l'operateur d'antisymetrisation A qui est etudie au paragraphe suivant. 

6.2.6 Action des operateurs S 2 et S z sur les determinants de Slater. 

L'action de l'operateur S z sur un determinant de Slater est de multiplier cd determinant par (n„ — 
np)/2, ou n a et rip designe les nombres de spinorbitales a et (5 

S z |ft 1 (l),A 2 (2), . . . , k N (N)\ = *,(1)|A 1 (1),A 2 (2), . . . , k N (N)\ + *,(2)|*i(l), * 2 (2), . . . , k N (N)\ + ... 

+«,(JV)|Ai(1),A 2 (2),---,*jv(^)| 
= - np)\k 1 (l),k 2 (2), k N (N)\ (6.33) 

Pour faire agir S 2 on utilise les relation 6.19, 6.20 et en general on obtient une combinaison lineaire 
de determinants de Slater. Si un determinant de Slater est toujours fonction propre de S z , il n'est pas 
toujours fonction propre de S 2 . Pour obtenir des fonctions propre de S 2 , on doit faire des combinaisons 
lineaires de determinants de Slater. Dans le cas des sytemes a couches completes, les orbitales (fonctions 
d'espace) sont toutes doublement occupees, c'est a dire qu'a chaque spinorbitale a, 4>j,, correspond une 
spinorbitale 4>j, de spin /? ayant la meme partie spatiale. 

D = |<M1)<M2) • ■ - 1)*jv(2JV)| (6.34) 

et 

S Z D = 

S 2 D = S+S-D = S-S+D = (6.35) 

en effet 

S+D = |^i(l)^i(2) . ..4>n{2N - 1)*jv(2JV)| + . . .+ |<M1)<M2) • -4n{2N - 1)* N {2N)\ = 
S-D = |^i(l)^i(2) . . .<j> N {2N - 1)$ N (2N)\ + . .. + 1^(1)^(2) . ..j N (2N - l)* N (2N)\ = (B.36) 

parceque les determinants qui apparaissent dans les sommes ont tous deux spinorbitales identiques, done 
deux lignes identiques et sont done nuls. Le determinant d'un systeme a couches completes est fonction 
propres de S 2 et de S z avec les valeurs propres S = 0, Ms = 0, c'est un etat singulet. 

On appele systeme a couches ouvertes, un systeme ou N orbitales sont monooccupees. Pour construire 
les fonctions propres de spin, il suffit de raisonner sur les 2^ determinants construits sur orbitales sim- 
plement occupees. Le determinant ne comprenant que des spinorbitales a est fonction propre de L z avec 
la valeur propre Ms = N/2, il est egalement fonction propre de L 2 puisque : 



S + \<f> 1 (lW)---<f>N(N)\ = 



(6.37) 



ce qui entrame 

s 2 \Mi)M2) ■ ■ -<Pn(n)\ = (s.s+ + s 2 + s z )\M^W) ■ ■ -4>n(n)\ = y(y + 1)1^(1)^(2) . ..<p N (N)\ 

(6.38) 

dS = 7 

Pour trouver les autres composantes du multiplet on fait agir l'operateur S- de facon successive. 
Les fonctions propres de S 2 correspondant au valeurs de S plus petites sont determiner en formant des 
combinaisons lineaires de determinants comme dans l'exemple suivant. 

Systeme a trois couches ouvertes. 

On peut former N 3 — 8 determinants : 
M s = | Di — |^i(l)^(2)^ 3 (3)| 

M s = \ D 2 = |^i(l)^(2)^ 3 (3)| £» 3 = |^i(l)^ 2 (2)^ 3 (3)| £» 4 = |^i(l)^ 2 (2)^ 3 (3)| 

M s = -\ £> B = |^i(l)^(2)^ 3 (3)| £» 6 = |^i(l)^ 2 (2)^ 3 (3)| I>7 = |tfi(l)^2(2)^ 3 (3)| 



M s = -\ D 8 = |^i(l)^2(2)^3(3)| 
Les determinants D\ et Ds sont fonctions propres de S 2 avec la valeur propre S = 



S 2 D 1 = (S-S+ + S 2 + S Z )D! = + |l>i + |^i = 

S 2 D 8 = (S+S- + S 2 - S 2 )D 8 = + ^D a + I As = ^-D 8 (6.39) 
En faisant agir S- on obtient : 

6(|-i) = ie(|,i) = -L(D 5 + D 6 + D 7 ) (6.40) 

Les composantes du doublet S = | sont calculees en formant des combinaisons lineaires de determinants 
orthogonales a ©(§, |) et a 0(|, — |) 

®{\,\) = aD 2 + bD 3 + cD A 
Q (\>-\) = dD 5 + eD 6 +fD 7 (6.41) 

comme 

s + e(\,\) = 

S -®{\'~\) = ( 6 - 42 ) 
= 

d + e + f = (6.43) 



Un choix possible est 



4'-^) = ^Pe-A) (6.44) 



o.z.i ^»peraieur a aniisymeirisaiion ei regies ae aiaier 



Le groupe des permutations S 



N- 



L'ensemble des n permutations de n objet forme im groupe appele groupe symetrique et note Sjv- Les 
permutations de deux objets sont appelees transpositions. Le groupe Sjv contient n(n — 1)/2 transpositions 
notees (a b). Les n— 1 transpositions (1 2), (1 3), . . . , (1 n) sont les generateurs du groupe. Une permutation 
paire est decomposable en un nombre pair de transpositions, une permutation impaire Test en un nombre 
impair. 

Exemple : le groupe S3 

Le groupe S3 contient 3! = 6 elements 



Pi 
P2 
P 3 
Pa 
P 5 
P 6 



= (1)(2)(3) 
= (1 2) 
= (1 3) 
= (2 3) 
= (12 3) 
= (13 2) 



(1 2)(1 2) 



(1 2)(1 3)(1 2) 
(1 2)(2 3) 
(1 3)(2 3) 



(1 3)(1 3) 



pair 

impair 

impair 

impair 

pair 

pair 



(6.45) 



L'operateur d'antisymetrisation. 

L'operateur antisymetriseur A est defini par la relation 



A=-^=Yt P P 



(6.46) 



la somme est effectuee sur les N permutations P du groupe Sjv, £p = +1 pour les permutations paires, 
ep = — 1 pour les permutations impaires. Dans le cas du groupe S3 on a : 



A = -j=[Pi -P2-P3-P4 + P5 + Pe] 



(6.47) 



Nous mentionnerons ici les proprietes de A qui seront utiles par la suite : 
1. Si Q est un element de Sjv alors : 



QA = -L= tpQP = CrR = e Q A 



(6.48) 



p 



2. 

A 2 - 

3. A est hermitique A^ — A 

Nous pouvons maintenant voir pourquoi A est appele antisymetriseur. Soit la fonction a TV electrons 
2, ...,N) obtenue par action de A sur $(1,2, . . . , N) 



Q*(1,2,...,JV) = QA*(1,2,...,N) 
= e q A*{l,2,...,N) 
= £Q *(1,2,...,JV) 



(6.50) 



itegies ae Slater 



Pour des determinants de Slater construits sur des spinorbitales orthogonales les valeurs des elements 
de matrice des operateurs monoelectroniques Q(i) et bielectroniques G(i,j) sont donnees par les regies 
de Slater. 

L'operateur H est la somme de contributions monoelectroniques (energie cinetique interaction cou- 
lombienne electron-noyaux) et bielectroniques (repulsion electronique). Tous les operateurs ayant une 
signification physique sont symetriques par rapport aux permutations des electrons. Si Q est un opera- 
teur quelconque ayant cette propriete il commute avec A: 



[-4,0] = [0,-4] = 
L'element de matrice de l'operateur Q a pour expression: 

(D K \Q\D L ) = J{Ak*(l)...k* N (N)}Q{Ah(l)...l N (N)}dr 

= J {k*(l) . ..k* N (N)}A^Q{Ah(l) . ..l N (N)}dr 

= J {k*(l) . ..k* N (N)}QA 2 {h(l) . . .l N (N)}dr 

= Vmj {k\{\) . ..k* N (N)}Q{Ah(l) . ..l N (N)}dr 

{*!(!) . . .k* N (N)}Qdet{h(l) . . J N (N)}dT 



(6.51) 



-h 



(6.52) 



Considerons d'abord le cas ou l'operateur Q est une constante que Ton pourra prendre egale a 1. L'element 
de matrice Qkl es t alors l'integrale de recouvrement entre deux determinants de Slater : 



Skl = J k* 1 (l)...k* N (N)det{l 1 (l)...l N (N)}dr 



(6.53) 



La multiplication de tous les elements d'une ligne par un facteur multiplie la valeur de ce determinant 
par ce meme facteur. Les facteurs &*(!) . . .k^{N) peuvent etre attribues a chaque ligne de sorte que : 



/ 



Skl = / DkiAt 



D 



KL 



*3(2)i!(2) *3(2)i 2 (2) 
k* N (N)h(N) k* N (N)l 2 (N) 



d'ou Ton tire 



Skl 



(h\h) (h\l 2 ) 
(k 2 \h) (k 2 h) 

(k N \h) (k N \h) 



*3(2)ijv(2) 
k* N (N)l N (N) 



(*iM 

(k 2 \l N ) 
(k N \l N ) 



(6.54) 



(6.55) 



(6.56) 



Dans cette expression (ki\lj) represente l'integrale de recouvrement des spinorbitales et lj. 

Pour calculer les elements de matrice d'operateurs monoelectroniques on utilisera le developpement 
de Laplace du determinant Dkl '■ 



N 



(6.57) 



m= 1 



oil les k- [i )l m [i ) sont les elements de la 1 ligne et de la m 
correspondants. -D^l est line fonction de toutes les variables electroniques a l'exception de celles de 
l'electron i. 

Si Q est im operateur monoelectronique de la forme : 



JV 



Q = E«W (6.58) 



t=i 

alors : 



JV JV „ JV JV 



(D K \Q\D L ) = E E / k*(i)g(i)l m (i)DT L dr = E E ( k i\lK)ST L (6.59) 



i — lm — 1 i — 1 m=l 



(ki\q\l m ) = J k*{l)q{l)l m {l)d Ti (6.60) 
En appliquant le developpement de Laplace aux cofacteurs D l ^ L de Dx l il vient : 

dtl= e ^Wnmfr (6-6i) 

n = l 



d'ou : 

JV 



k*{i)i m {i) k*{i)i n {i) 

k*(j)l m (j) k*(j)l n (j) 



N 

= e *;(o*;(j){i m (0' n (j) - ^urn^r m 



Si Ton considere l'operateur bielectronique symetrique : 

JV 



G = ^j(i,j) (6.63) 



on obtient : 

JV JV 



Pjc|gii> l > = EE /*. T (o*;o')i7(i,j){' m (0'n(j)-'n(0' m (j)}^r^ 

JV JV 

= ^^[(MjMWnl-^ijlslUmKr ( 6 - 64 ) 

ou : 

(MjM^rJn) = [WmM&^n] = ^ k* (l)k* g(i, j)l m (i)l n (j)dTidTj (6.65) 

Les expressions ci dessus se simplifient considerablement si les determinants Dk et sont construits 
sur une base orthonormalisee de spinorbitales. Dans ce cas les integrales de recouvrement (&,• \lj) sont nulles 
si les spinorbitales et lj sont differentes. On peut ordonner le determinant de sorte que toutes les 
spinorbitales communes avec Dk soient situees dans les memes colonnes dans les deux determinants, ce 
qui entraine : 

(*i|ij) = Vi^i (6.66) 

Le determinant des recouvrements Skl es t alors diagonal; ses elements diagonaux sont egaux a 1 pour 
les spinorbitales qui coincident et nuls pour celles qui ne coincident pas. Les determinants construits sur 



aes spmorDitales ortnonormees sont normes, s lis sont aistmcts Us sont ortliogonaux entre eux. Lie plus 
tous les cofacteurs non diagonaux sont nuls et : 

= AA m (6.67) 
avec A,- = 1 si Dk et ne different que par leur i eme spinorbitale, sinon Aj = 0. De meme : 

Sk'l" - AijSimSjn (6.68) 

avec A,j = 1 si Dk et ne different que par leur i eme et j eme spinorbitale, sinon A,j = 0. 
L 'expression des elements de matrice mono et bi-electroniques deviennent respectivement : 

N 

{D K \Q\D L ) = ^2{k i \q\l i )A i (6.69) 
t=i 

et 

N 

(D K \G\D L ) = J^Kkikjlglhlj) - (kikjWjWAij (6.70) 

i<3 

Pour que les elements de matrice ne soient pas nuls il y a trois possibilites : 

1. D K EE D L 

N 

{D K \Q\D L ) = '£ l {k i \q\l i ) (6.71) 

t=i 

N 

(D K \G\D L ) = YftWMiih) ~ {kikjlglljh)} (6.72) 

i<3 

2. Dx et Dl different par la i eme spinorbitale 

(D K \Q\D L ) = {ki\ q \k) (6-73) 

N 

(D K \G\D L ) = Y^K^Mhh) ~ {kikjlglljh)} (6.74) 

3. Dk et Dl different par deux spinorbitales 

(D K \G\D L ) = {kikjlglhlj) - {kikjlglljh) (6.75) 



Chapitre 7 

La molecule de dihydrogene H2 



L'hamiltonien electronique de la molecule de dihydrogene a pour expression : 

g = -Iy;-Iy» 1 --- 1 --- 1 ,-, 1 , +r JL_ + , 1 , (7.1) 

2 1 2 2 |ri-R^| |ri-R B | l^-R^I |r 2 -R B | |r 2 — n | |R B -R 4 | V ' 

7.1 La fonction d'onde de Heitler et London 

La molecule de dihydrogene possedant deux electrons on peut ecrire sa fonction d'onde sous la forme 
du produit d'une fonction d'espace et d'une fonction de spin construite sur les fonctions a et /?. 

¥(n , r 2 , <r u <r 2 ) = *(n, r 2 )0(S, M s ) (7.2) 

Une solution approchee simple est la fonction de Heitler et London : 

, , ls^(ri)ls B (r 2 ) zfc ls^(r 2 )ls B (ri) 

V ' V2(l ± S 2 ) V ^ 

ou Isa et Isg designent des fonctions Is centrees sur les noyaux A et B et S l'integrale de recouvrement : 

S = J \ SA {v)\ SB {v)dv (7.4) 

Si Ton permute les coordonnees d'espace des deux electrons 

$ + (r 2 ,ri) = $ + (ri,r 2 ) symetrique 

$ + (r 2 ,ri) = — $ + (ri,r 2 ) antisymetrique (7-5) 

Les fonctions totales satisfaisant le principe d'antisymetrie sont d'une part le singulet 

% = ^+ M*i)/?M - (7.6) 
et les trois composantes du triplet 

^(a^O^J + a^a)^)) (7.7) 

/%l)/%2) 

La fonction de Heitler et London tend vers la fonction d'onde exacte quand la distance internucleaire R 
tend vers l'infini. 
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i.i.i evaluation ae i energie 

L'energie du systeme a pour expression 



E ± = {$ ± \H\$ ± ) = -l + f^ 



n - R? 



1 1 1 
+ i r + 



Q = J J l^( ri )l SB (r 2 ) 

f f r i i i i 

A = / ls A ( ri )l SB (r 2 ) -- ==-,-1- ^n + r- ri+M 

J J I |ri — Kb| \r 2 — H A \ |r2 — ri| IRb ~ 



r 2 - R.a I |r 2 — rij \R B - R A \ 



(7.8) 

ls,4 (ri)lsB(r 2 )dridr 2 

ls B {ti)\s a (r 2 )dridr 2 
(7.9) 



ou en detaillant 



Q = - [ U *te\ dT 2 - [ -MlgL dri + / / U A ( Tl )U A ( Tl ) 1 , l 8B (r 2 )l 8B (r 2 )dndr 2 
J |r 2 — K4I J |ri — K B | J J |r 2 - ri| 



+ ■ 



R« — R, 



J |r 2 -R J i| J |ri-R_e| 

+ / / ls,i(ri)lsB(ri) 1 f lsyi(r 2 )ls B (r 2 )dridr 2 + — — r 

J J |r 2 - ri| |R B - R^| 

L'integrale 

(aa\bb) = / / ^^(r^ls^^i), j4s B (r 2 )lsB(r 2 )dridr 2 

J J l r 2 - ri| 

est appelee integrale coulombienne, et 

(ab\ab) = f f 1s a (ti)1s b (ti)- f lsyi(r 2 )ls B (r 2 )dridr 2 

J J l r 2 - ri| 



(7.10) 
(7.11) 

(7.12) 



integrale d'echange. Dans la fonction d'onde de Heitler-London, l'exposant des orbitales ls^ et ls B 
est egal a 1.0, si l'exposant des fonctions de base est pris comme parametre variationnel, on ameliore 
sensiblement les resultats, et en particulier la distance internucleaire d'equilibre, comme le montre la 
tableau 7.1 





C = i.o 


C = 1-166 


exact 


D e (e.V.) 


3.20 


3.78 


4.75 


R e (u.a.) 


1.64 


1.41 


1.40 



Tab. 7.1: Energie de dissociation et distance d'equilibre calculees par les fonctions d'Heitler-London, sans 
et avec optimisation de l'exposant 

La fonction d'Heitler et London ne peut pas etre ecrite sous la forme d'un determinant de Slater, elle 
appartient a un type de fonctions approchees appelees valence-bond (VB). 

Les fonctions VB $ + et font apparaitre des produits de fonctions de base, centrees sur des noyaux 
differents et sont pour cette raison appelees fonctions covalentes. En plus de ces deux fonctions covalentes 
ont peut former deux fonctions ioniques : 



$i±(ri,r 2 ) 



1s a (ti)1s a (t 2 ) ± lsB(r 2 )ls B (ri) 



(7.13) 



qui correspondent aux situations ou les deux electrons sont localisees sur 1 un des noyeau. 

La fonction covalente ne possede pas de plan nodal, elle est de symetrie et correspond a l'etat 
singulet 1 S+ la fonction possede un plan nodal et correspond au triplet 3 £j. 



7.2 Methode des orbitales moleculaires 



Pour les orbitales moleculaires les plus simples construites a partir de fonctions de base centrees 
sur les noyaux sont : 



<Pg( r ) = 
fu(r) = 



U a (t) + U b (t) 
ls A (r) - ls B (r) 



(7.14) 



A partir de ces orbitales on construit les spinorbitales 



et les determinants de Slater 

Di = 

D 2 = 

D 3 = 

D A = 

D 5 = 

D 6 = 



1 
1 
1 

1 
1 
1 

~7i 



-M 1 ) = <Pg( r i) a (<ri) 
4>„(1) = <p u (ri)a(ai) 

^(1) = p g (riM<ri) 



1 



(7.15) 



4> B 


-)<M2)I 


4>g0 


0^(2)1 


4>g0 


0^(2)1 


tgO 


0^(2)1 


tgO 


-)<M2)I 


<t>u{ 


0^(2)1 



(7.16) 



a partir desquels on forme les fonctions d'onde de trois etats singulets 



^1(1, 2) = D 1 
^2(1,2) = -L(D 3 -D 4 ) 
^3(1,2) = D 6 
et les trois composantes du triplet 3 £+ 

3 *i(l,2) 
3 * 2 (1,2) 

3 * 3 (1,2) 



etat fondamental singulet 1 Sj" 
etat excite singulet 1 !] J 
etat excite singulet 1 !]^ 



(7.17) 



D 2 
1 

D 5 



(D 3 + D 4 ) 



(7.1f 



i.a.l evaluation ae i energie 

L'operateur hamiltonien du systeme peut etre decompose en quatre termes 



H = hS\) + h 1 (2) + — + 4 (7.19) 



h\i) = -lyi 



2 |ri-R 4 | |ri-R B | 

ft '< 2 > = (7 - 20) 

Les energies des differents etats sont donnees par les relations : 

£iCS+) = i ((£> 3 |^|£> 3 >-2(£> 3 |^|I>4> + P4|^|£>4> 



^CS+J = (£> 6 |^|£>, 



^s( 3 S+) = ^ (( J D 3 | J ff| J D 3 } + 2{D 3 \H\D 4 ) + {D A \H\D A )) (7.21) 
et evalues en fonction des integrales sur les spinorbitales a l'aide des regies de Slater : 

SoCS+J = (^l/i 1 !^) + (^lA 1 !^) + MaMg) " M a \4> a 4>g) = 2(^1^1^) + MaMg) 

^CS+J = 2(^|/i 1 |^) + (^^|^^) 

£ 3 ( 3 £+) = ^ B |A 1 |^ B > + ^ U |A 1 |^> + (^ U |^ U > + (^ U K^> (7-22) 
qui sont a leur tour evaluees sur les fonctions de base ls^ = a, \sb = b, avec cette convention : 

(^l^) = ^[H/^W + H/^)] 

(tulh 1 ^) = T ^[(a|/i 1 |a)-(a|/i 1 |6)] 

{<f>g<f>g\<f>g<f>g} = on \ Q -, 2 [(aa\aa) + 4(aa\ab) + (aa|66) + 2(ab\ab)] 
2(1 -\- b) 



[(aa\aa) + (aa\bb) — 2(ab\ab)] 



2(1 - S 2 ) 



{<f>g<f>u\<f>u<f>g} = 2(1 - S 2 ) [(""I"") ~ ( aa \ bb )] 

{<f>u<f>u\<f>u<f>u} = 2(1 - S) 2 ^ aa l aa ) ~ 4 ( aa l afc ) + (aa\bb) + 2(ab\ab)] (7.23) 

ou 

(alh^b) = J afriJh^nWrOdri 

(a6|cd) = j j a'(ri)6(n)^-c'(r 2 )d(r 2 )dridr 2 (7.24) 

La notation (a6|cd) des integrales bielectroniques fait apparaitre l'interaction de deux distributions de 
charges, une autre notation est egalement utilisee : 

[ab\cd] = j ya'(ri)6'(r2)^-c(ri)d(r 2 )dridr2 (7.25) 



Le grapne de 1 energie de ces etats en tonction de la distance mternucleaire It est presente ngure (.1 La 
courbe representant l'etat fondamental 1 S+ presente son minimum pour R = 1.6 u.a., valeur voisine de 
la valeur experimentale 1.4 u.a. Cependant l'energie de dissociation calculee D e = 2.65 eV ne represents 
que 55% de la valeur experimentale, et est inferieure a la valeur calculee avec la fonction de Heitler et 
London. Aux grandes valeur de R, l'erreur sur l'energie est bien plus grande (8.5 eV pour R = oo) et le 
calcul predit la dissociation ionique H + + H~. Les courbes des etats 1 Sj et 3 Ej sont qualitativement 
correctes. L'energie d'excitation calculee pour l'etat singulet, T e = 13.0 eV est en accord satisfaisant avec 
sa valeur experimentale 11.4 eV. 

On peut ameliorer de facon significative la fonction d'onde approchee de l'etat fondamental en op- 
timisant l'exposant des fonctions Is. Avec (* = 1.197, l'energie de dissociation est D e = 3.488 eV et la 
distance d'equilibre R = 1.38 u.a. 




Fig. 7.1: Courbes d'energie potentielle des etats 1 ^ g , 1 et 3 S U de la molecule H2 



f.z.z uomparaison aes ioncnons vd ei ivivj 

Les fonctions MO decrites preceedement peuvent etre exprimees sous la forme de combinaisons li- 
neaires de fonctions valence-bond. Ainsi : 

i(l, 2) <x (1^(1) + U B {1)) (1^(2) + 1« B (2)) (a(l)/?(2) - /?(l)a(2)) 

<x [(1^(1)1^(2) + U B (1)U B (2)) + (l^(l)l« fl (2) + l« fl (l)l^(2))] (a(l)/?(2) - /?(l)a(2)) 
<x [<D i+ + <D + ] (a(l)/?(2) - /?(l)a(2)) (7.26) 

La fonction d'onde MO contient une contribution de la fonction valence bond ionique, ce qui explique 
pourquoi l'energie de dissociation calculee par la methode des orbitales moleculaire est inferieure a la 
valeur valence bond. 
De meme : 

^2(1,2) <x [* 4 _ -*_](a(l)/?(2)-/?(l)a(2)) 

1 M/ 3 (1,2) <x [$,■+ - <&+] (a(l)/?(2) - /?(l)a(2)) 

3 *i(l,2) oc $_a(l)a(2) (7.27) 



7.3 Interaction de configuration 



Au lieu de construire la fonction d'onde sur un seul determinant ou sur une seule fonction VB on peut 
utiliser une combinaison lineaire de fonctions de meme symetrie. Les coefficients du developpement sont 
de nouveaux parametres variationnels. Ainsi, avec les fonctions d'espace VB on exprimera la fonction 
d'espace IC de l'etat fondamental par la combinaison lineaire : 

$ = Cl $ + + c 2 $ i+ (7.28) 

Pour determiner les coefficients c\ et c 2 on doit minimiser la fonctionnelle 



En posant 



cjSn + 2cic 2 Si2 + c\S 22 



soit 



(7.29) 



H n = <$+|ff|$+> 

H 12 = H 21 = (<S> + \H\<S> i+ ) 

H 22 = {$ i+ \H\$ i+ ) 

Sn = ($+|$+) 

Si 2 = S 2 l = ($ + |$ i+ ) 

S22 = (7.30) 



(7.31) 



E (cfSn + 2c lC2 Si2 + c\S 22 ) = c\R lx + 2 Cl c 2 H 12 + c\H 22 (7.32) 
Les valeurs optimales des coefficients correspondent a la condition variationnelle 



en aerivant par rapport a c\ et c 2 

B E 

2E(c 1 S 11 + c 2 S 12 ) + (c 2 S 11 +2c 1 c 2 S 12 + c 2 2 S 22 ) — = 2ci#n + 2c 2 H 12 

f) E 

2E(c 1 S 12 + c 2 S 22 ) + (c 2 S 11 +2c 1 c 2 S 12 + c 2 S 22 ) — = 2 Cl H 12 + 2c 2 H 22 (7.34) 

Cc 2 

soit, en tenant compte de la condition variatinnelle 7.33 

(H n - ES^d + iH^- ES 12 )c 2 = 

{H 12 - ES 12 ) Cl + (H 22 - ES 22 )c 2 = (7.35) 

En apparence, ce systeme est indetermine parcequ'il y a trois inconnues, E, c\ et c 2 pour deux equations. 
En fait c\ et c 2 ne sont pas independents en raison de la normalisation 

c\Sn + 2c lC2 S 12 + c 2 S 22 = 1 (7.36) 

Pour resoudre on prend le rapport c 2 /c\ comme variable auxilliaire 

ci _ Hn — ESu _ Hi 2 — ESi 2 ^ ^ 

c 2 H\ 2 — ES\ 2 H 22 — ES 22 

ce qui impose 

(iTn- ESn)(H 22 - ES 22 )-(H 12 - ES 12 f = (7.38) 
La resolution de cette equation est equivalente a la resolution du determinant seculaire 



#n — -SSii #i 2 — ES\ 2 
H\ 2 — ES\ 2 H 22 — ES\ 



(7.39) 
ce qui 



En utilisant des fonctions VB a exposant optimises, on trouve R e = 1.43 u.a. et D e = 4.025 eV, 
constitue une amelioration des resultats obtenues avec une fonction VB simple. 

7.4 La fonction d'onde de James et Coolidge 

Les differentes fonctions d'onde decrites precedement divergent toutes pour ri 2 = 0, la nature de cette 
singularite peut etre analysee en ne considerant que les termes important de 1'hamiltonien quand ri 2 — > 0, 
ce qui correspond a l'equation de Schrodinger approchee 

lv i - l-yl + — ) $ = E ® ( 7 - 40 ) 



2 1 2 z n 2 

que Ton peut transformer en introduisant les coordonnees du "centre de gravite" 

R = ^ (ri + r 2 ) r = ^ (ri - r 2 ) r = r 12 (7.41) 



et apres separation des mouvements 



1 1 > 

+ - ) $ = #$ (7.42) 



2 



En developpant $ en serie de r 

$ = a + br + cr 2 + . . . (7.43) 
En substituant dans l'equation 7.42 et en s'arretant au premier odre on trouve 

b = l -a (7.44) 



Une ionction a onae approcnee de tres Donne qualite doit done contenir des termes imeaires en r\i. James 
et Coolidge ont optimise une fonction en coordonnees elliptiques de cette forme 

$(1,2) = i- exp (-0.75(6 + 6)) {2.23779 + 0.80483 (ij? + iji) 

-0.559966 - 0.60985 (6 + 6) + 0.56906ri 2 } (7.45) 

Cette fonction permet de calculer une energie de dissociation D e = 4.53 eV. La valeur experimentale non- 
relativiste a ete calculee par Kolos et Wolniewicz avec une fonction du type James Coolidge comprenant 
80 termes. 



Chapitre 8 

La methode Hartree-Fock 



8.1 Approximations de champ moyen (electrons independants) 

8.1.1 Developpement de la fonction d'onde sur une base de determinants 

Une propriete importante des determinants de Slater est le theoreme d'expansion : les determinants 
de Slater d'ordre N construits sur une base complete de spinorbitales forment une base complete pour 
les fonctions antisymetriques a iV-fermions. Cette base de spinorbitales est formee par l'ensemble des 
fonctions propres d'un operateur monoelectronique qui reste a definir. Le probleme a resoudre est done 
de trouver un tel operateur. 

8.1.2 Invariance des determinants, non-invariance des spinorbitales 

Considerons deux ensembles de TV spinorbitales {tp} et {<£>'}, tels que les ip 1 soit des combinaisons 
lineaires des ip : 

<p' = <M (8.1) 



D= ^deUp 

D '= vk det< ^' 

En substituant tp par son expression en fonction de ip : 

1 



(8.2) 



D 1 = 



1 



= D x det A 



det ipA 
det <^?det A 



(8.3) 



Les deux determinants sont egaux a une constante multiplicative pres. lis representent la meme situation 
physique, de plus si ip et ip' sont tout deux normes, alors D et D' sont identiques. 



8.1.3 Equations Hartree-Fock et formulation LCAO de Roothaan 

Comme nous l'avons deja souligne, 1'hamiltonien electronique est independant des variables de spin 
electronique. La fonction d'onde, elle, depend a la fois des variables d'espace et des variables de spin. 
Cela a deux consequences : 



71 



1. \^) est tonction propre des operateurs ae spin total b et d z 



S 2 \ty = S(S + 1)\*S) 
S 2 \V) = M S \V) 

2. les electrons etant des fermions, la fonction d'onde doit etre antisymetrique par rapport a la per- 
mutation de deux particules quelconques : 

i , «|*> = -|*> 

Compte tenu de ces remarques la structure mathematique la plus simple pour une solution approchee 
de l'equation d'onde electronique est celle d'une somme de determinants de Slater construits sur des 
spinorbitales. 

\^}^J2 CkDk ( 8 - 4 ) 

K 

les coefficients Ck sont determines a partir des contraintes liees au spin et a la symetrie d'espace, pour 
un systeme a N electrons les determinants de Slater Dk ont pour expression : 



D 



K 



4> kl {\) ^ 2 (i) ... 4> k „{\) 

4> kl {2) <^ 2 (2) ... ^(2) 
MN) 4>k a '(N) ... <f> kN (N) 



(8.5) 



Exemple : l'etat fondamental de H2 



1 
1 



•Mi) 
<M2) 



^i(i) 
<M2) 



- 7 ={^ 1 (l)^ 1 (2) 



tfi(2)^(l)} 



(8.6) 



En raison des contraintes de symetrie d'espace les spinorbitales ont la meme partie spatiale. 

La derivation des equations Hartree-Fock peut d'abord etre effectuee dans le cas d'un seul determinant 
construit sur des spinorbitales orthonormees, puis generalisee au cas des combinaisons lineaires de deter- 
minants. La derivation esquissee ici correspond a ce premier cas (UHF). 

Considerons une fonction d'onde approchee de la forme : 

l*> « \D) 

le principe variationnel stipule que ^jjjpy^ est une borne superieur de l'energie de l'etat le plus stable 
de la symetrie considered : 

(D\H\D) 



(D\D) 



> E 



La determination des spinorbitales peut done etre effectuee en minimisant (D\H\D) tout en maintenant 
les spinorbitales orthonormees. On utilisera done une methode de minimisation avec contrainte ou sont 
introduits des multiplicateurs de Lagrange. La fonctionnelle a minimiser est : 



N 



{D\H\D)-J2^m 



(8.7) 



L application cm prmcipe vanationnel a cette ionctionnelle s ecnt 



N 



S\{D\H\D)-Y,^m\ (8- 



N 



8{D\H\D) = ^{(^l^l^ + ^lft 1 ^}} 

i = l 

N N r i i 

+ ] ( 5 fifi\ — \fifi) - { s <Pi<Pi\ — \<Pi<Pi) 

\^-\fj S( Pi)\ 
ri2 ) 

Y J {^i\F\fi) + {fi\F\^i) (8- 



i=i j=i 



+ {fifi\ — \ s fifi) ~ {<Pi<Pj\ 
ri2 

N 



i = l 



On montre facilement que 



JV 



S(D\H\D) = J2(^i\F\<f>i) + (MHHi) (8.10) 



i = l 

ou i* 1 est un operateur monoelectronique : 

^=-^-E7t+i:^- A 'i 
a ^ i=i 
- — „ — ' 

/.MO 

qui contient la contribution cinetique et potentielle electron-noyaux et ou la repulsion electron-electron 
est remplacee par sa moyenne. Les operateurs de Coulomb, J , et d'echange, K , ont pour expression : 

Jj(lMl) = (| ^(2)-^^(2)dr 2 ) tf(l) (8.12) 

tfj (1)0(1) = (| <^(2)-^(2)dr 2 ) ^(1) (8.13) 
Finalement on obtient un systeme de TV pseudo-equations aux valeurs propres : 

F\4>i) = *Mi) (i = l,N) (8.14) 

L'operateur F est un operateur hermitique, ses fonctions propres torment une base complete et infinie. 
l'energie totale du systeme est donnee par : 

E?I = ^{ilh'li) + ljr{{ij\ij) - + £ ^ (8.15) 

t = l j = l >1<B AB 

Ces equations sont facilement generalisables aux cas (RHF) ou les contraintes de spin et de symetrie spa- 
tiale sont prises en compte. En particulier pour les systemes a couches fermees, le cas de l'etat fondamental 
de la plupart des molecules, la contrainte de symetrie impose qu'a chaque spinorbitale a corresponde une 
spinorbitale /3 ayant la meme partie spatiale. Le determinant s'ecrit alors : 

1 TV 
D = —7= det{^i^i^ 2 <^2 • • • <Pm<Pm) avec M = — (8.16) 

viV ! ^ 



1 Uperateur de fock aevient : 

M 

F = h 1 +Y J 2j i- k i ( 8 - 17 ) 

et Ton doit resoudre M equations de la forme : 

F\4>i) = uHi) (8.18) 

Dans le cas des systemes a couches ouvertes les equations RHF sont plus difficiles a etablir. Souvent on 
se contentera d'une approximation UHF sans contrainte de spin. 

8.1.4 Resolution des equations Hartree-Fock 

L'operateur de Fock depend de ses solutions au travers des operateurs coulombien et d'echange. Les 
equations Hartree-Fock forment un systeme d'equations integro-differentielles couplees qui doivent etre 
resolues simultanement. Dans la pratique on utilise une methode iterative appelee methode du champs 
autocoherent, en ang lais SCF (Self Consistent Field). A l'iteration l'operateur est construit a partir 
d'un jeu d'orbitales d'essai |<^ )}. Ce premier fockien est utilise pour calculer un nouveau jeu d'orbitales, 
|<?)W} qui sera utilise pour construire F^'l et ainsi de suite jusqu'a l'obtention d'une bonne convergence 
sur l'energie (les valeurs propres, la matrice densite) : 

pin-i)^) = e\ n W>) (8.19) 

Dans la pratique les equations Hartree-Fock sont resolues de facon approchee en developpant les solutions 
sur une base finie, done incomplete, de fonctions connues {x}- Ces fonctions sont en general centrees sur 
les atomes et sont pour cette raison appelees orbitales atomiques de base. Sous forme matricielle on a : 

* = xc 

F X C = e X C (8.20) 



En multipliant a gauche par '■ 



En multipliant a gauche par '■ 



et 



il vient : 



X^F X C = ex^xC 

F X C = e X C (8.21) 

X 1 F X C = ex^xC (8.22) 

X 1 X = S (8.23) 

X 1 F X = F (8.24) 



FC = eSC (8.25) 
oil e est la matrice des energies orbitalaires. En posant : 

D = S 1 / 2 C ou C = S" 1 '' 2 D (8.26) 



on trouve 



s -i/2 FS -i/2 D _ £ g- 1 / 2 S 1/2 D = £D 



(8.27) 



Un est alors ramene a la aiagonalisation a line matrice symetrique. 

Pour les systeme a couches fermees, l'expression explicite des elements de la matrice de Fock en 
fonction des integrates sur les fonctions de base est donnee par la relation : 



X.o 



{liv\\a) - -{li\\va) 



(8.28) 



(H Acr ) = J J x^(ri)#V4 r i)^XA(r2)#V<7(r 2 )dridr 2 

occ 



(8.29) 



Les Pp,, sont les elements de la matrice densite. L'energie totale a pour expression : 



H,V A,<7 



(HAer) - ^(pM^) 



A<B 



Rab 



.30) 



8.1.5 Limitations de la methode Hartree-Fock 

1. La methode Hartree-Fock est une methode variationnelle. Elle permet de calculer une fonction 
d'onde approchee pour l'etat fondamental mais un calcul de meme qualite est impossible pour les 
etats excites de meme symetrie et de meme multiplicite de spin. 

2. La methode Hartree-Fock favorise la dissociation ionique des molecules. Cela a deux consequences : 
les distances internucleaires entre atomes lies sont calculees trop courtes, les constantes de force 
d'elongation sont surestimees d'environ 20%. ceci peut etre montre sur le cas de H 2 

Les operateur de Fock pour H 2 , H~ et l'hamiltonien de H sont respectivement : 



F = 




1 


1 






TiA 


TiB 


F = 




1 


+ J 


H = 




1 

TiA 





(8.31) 



Pour une distance Rab tendant vers l'infini, l'operateur de Fock de H 2 tend vers celui de H + H + 
et non vers l'hamiltonien de deux atomes d'hydrogene sans interaction. 

3. La methode Hartree-Fock est une methode a particules independantes ce qui veut dire que deux 
electrons de spins antiparalleles peuvent occuper simultanement la meme position. Cela a pour 
consequence d'augmenter la repulsion coulombienne electron-electron de maniere non physique. 
Une fonction d'onde satisfaisant les conditions aux limites physiques doit s'annuler lorsque plusieurs 
particules occupent la meme position. II existe autour de chaque particule un volume dans lequel la 
probabilite de trouver une autre particule est negligeable. Deux causes physiques peuvent inter venir : 

(a) la repulsion coulombienne qui explique la formation du trou de Coulomb 

(b) la repulsion de Pauli entre particules de spin parallele responsable de la formation du trou de 
Fermi qui est pris en compte par l'approximation Hartree-Fock. 




Fig. 8.1: Energie de l'etat fondamental de H 2 



8.1.6 Qualite des calculs SCF 

On peut obtenir des fonctions Hartree-Fock exactes, on parle alors de limite Hartree-Fock, par des 
methodes numeriques limitees essentiellement aux atomes et aux molecules diatomiques. Tous les calculs 
effectues avec des bases finies donnent des energies superieures a la limite Hartree-Fock. C'est le choix 
de la base qui permet de presumer de la qualite d'un calcul SCF. Le choix de la base est l'etape cruciale 
dont depend le cout du calcul. Le temps de calcul est proportionnel a N£, Nf, designant la dimension 
de la base. Le cout des etapes post SCF est proportionnel a N^. II y a une limitation a la fois physique 
et economique a l'augmentation de la taille des bases. L'experience a montre que des fonctions de base 
optimisees sur les atomes isoles conviennent tres bien aux calculs molecul aires. On distingue deux types 
principaux d'orbitales atomiques : 

1 Les orbitales de Slater ou STO 

Xmm{r; C) = N^r"- 1 exp (-Cr)Y, m (^ 4>) 

Ces orbitales ont un comportement correct a la fois pour r — > et r — > oo. 

Selon le nombre d'orbitales de Slater utilisees pour decrire un niveau occupe de l'atome on distinguera 
les bases : 

- simple-^ ou minimales dont les exposants sont determines par les regies d'ecrantage de Slater ou de 
Clementi et Raimondi 

- double-^ comme les bases de Clementi et Roetti 

- triple-^ 

Les orbitales de Slater ont ete longtemps utilisees pour des calculs sur des molecules diatomiques ou 
lineaires (programme ALCHEMY). Le calcul des integrales multicentriques qui interviennent dans la 
construction de la matrice de Fock est tres laborieux et pour cette raison il n'existe pratiquement aucun 
logiciel les utilisant pour des calculs moleculaires. 

2 Les orbitales gaussiennes ou GTO 

Ximn (r; a) = N lmn (x - X A )'(y - Y A ) m (z - Z A ) n exp (-ar 2 ) 



permettent im calcul emcace aes mtegrales multicentriques et c est la raison prmcipale ae leur succes. 
Cependant leur comportement au voisinage du noyau aussi bien qu'a grande distance est incorrect. Pour 
atteindre une qualite semblable a celle d'un calcul effectue avec des orbitales de Slater, il faut une base 
beaucoup plus grande. Afin de limiter le nombre d'integrales on utilise des contractions de gaussiennes : 

n c 

Vlm„ (r) = ^2 CiXlmn {r\ <*i) 
t=l 

La terminologie utilisee pour les orbitales de Slater a ete etendue aux gaussiennes, mais dans ce cas 
"zeta" signifie contraction. Dans le cas des bases minimales, Pople a introduit la nomenclature STO-nG 
qui signifie que chaque contraction representant une orbitale de Slater est constitute de n primitives. 

On appelle base split- valence, une base ou les orbitales atomiques des cceurs sont decrites par une 
seule contraction, alors que les orbitales de valence le sont par au moins deux. 

On peut ameliorer la base en introduisant des orbitales de polarisation. Les orbitales de polarisation 
sont des orbitales dont le nombre quantique angulaire I est plus grand que ceux des orbitales occupees 
de l'atome, par exemple des orbitales p sur l'hydrogene, d sur le carbone. 

Criteres pour le choix des bases 

Le choix de la base depend essentiellement de ce que Ton desire calculer. On peut cependant enoncer 
quelques regies general es. 

1. Examen de la litterature sur la meme molecule ou sur des composes voisins. Theochem publie 
chaque annee une bibliographie complete des calculs ab initio. 

2. Homogeneite de la base. On doit veiller a ce que la qualite des fonctions de bases soit uniforme sur 
tous les atomes. 

3. Balance : si Ton introduit des fonctions de polarisation il faut en introduire sur tous les atomes. 

Quelques indications sur les bases les plus utilisees 

Les bases les plus utilisees sont les bases de Pople et al dans lesquelles les exposants des fonctions s 
et p d'une meme couche sont egaux. la nomenclature est la suivante : 
Bases minimales : STO-nG n — 3, 4, . . . , 6 

Bases split-valence : 

- n-21 G n = 3, 4,..., 6 

- n-31 G 

- n-311 G 

Les orbitales de polarisation sur les atomes lourds sont indiquees par une etoile, sur les hydrogenes par 
une seconde etoile : 
Exemple 6-31 G*, 4-31G** 

L'addition d'orbitales diffuses est signalee par un + : 6-31G**+ Les bases minimales STO-nG, prin- 
cipalement la base STO-3G, sont surtout utilisees au calcul de systemes comprenant plusieurs dizaines 
d'atomes, et pourlesquels on recherche des resultats qualitatifs concernant la structure. Les surfaces 
d'energie potentielles et les proprietes electriques calculees avec ces bases sont peu fiables. 

Les bases split-valence sont, aujourd'hui les bases les plus utilisees, pour la plupart des problemes 
de physico-chimie la presence d'orbitales de polarisation est necessaire des que Ton desire des resultats 
semi-quantitatifs. Les bases n-311 G sont plus specifiquement adaptees aux calculs tenant compte de la 
correlation electronique. Pour calculer des anions il est necessaire d 'introduire des orbitales diffuses 

Bases de Huzinaga et contractions de Dunning : (9s5p) / [4s3p], (10s6p) / [5s4p] 

Les bases MINI et MIDI d'Huzinaga. 

Les bases de van Duijneveldt (13s8p) 



uase 


^xi4 


imxi 3 


X12VJ 


in 


ST0-3G 


-39.727 


-55.454 


-74.963 


-98.571 


4-31G 


-40.140 


-56.102 


-75.907 


-99.887 


6-31G* 


-40.195 


-56.184 


-76.011 


-100.003 


6-31G** 


-40.202 


-56.195 


-76.023 


-100.011 


limite HF 


-40.225 


-56.225 


-76.065 


-100.071 



Tab. 8.1: Energies de molecules a dix electrons (a.u.) 



base 


CH 4 


NH 3 


H 2 


HF 


STO-3G 


2.047 


1.952 


1.871 


1.807 


4-3 1G 


2.043 


1.873 


1.797 


1.742 


6-31G* 


2.048 


1.897 


1.791 


1.722 


6-31G** 


2.048 


1.897 


1.782 


1.703 


limite HF 


2.048 


1.890 


1.776 


1.696 


Exp. 


2.050 


1.912 


1.809 


1.733 



Tab. 8.2: Distances d'equilibre (u.a.) 



Les bases de van Duijneveld sont un bon point de depart lorsque Ton veut approcher la limite Hartree- 
Fock. Dans ce cas il faudra introduire plusieurs orbitales de polarisation. Les performances de bases de 
differente qualite peuvent etre evaluees en examinant les tableaux 8.1, 8.2, 8.3. 

Les resultats reportes dans le tableau 8.1 montrent que l'energie converge assez rapidement vers la 
limite Hartree-Fock, le fait d'utiliser une base split-valence abaisse l'energie de 0.4 a 1.3 u.a., le gain 
obtenu en ajoutant des orbitales de polarisation est plus modeste, de l'ordre de 0.1 hartree. 

L'erreur sur les distances interatomique est toujours tres faible, inferieure a 0.1 A, avec des bases 
split- valence polarisees elle est de l'ordre du centieme d'angstrom. 

La valeur calculee du moment dipolaire est tres sensible a la nature de la base. En general, les bases 
minimales ont tendance a sous estimer le moment dipolaire, tandis que les bases split- valence le surestime. 
L'addition d'orbitale de polarisation a pour effet d'abaisser la valeur calculee du moment dipolaire et de 
la faire tendre vers la limite Hartree-Fock. Pour obtenir un accord excellent entre valeurs experimentales 
et calculees il faut faire intervenir la correlation electronique. 

II existe un artefact du aux bases incompletes lors du calcul de molecules en interaction : l'erreur de 
superposition ou BSSE. L'energie d'interaction d'un complexe AB est : 

-^int = Eab — Ea — Eg 

avec une base complete il n'y aurait aucun probleme mais les bases utilisees pour calculer Eab,Ea,Eb 
sont de dimension differentes done incompletes de facon differente. Dans le complexe, les molecules A et 
B sont mieux decrites que dans les monomeres parce que les orbitales de la molecule partenaire jouent le 
role d'orbitales de polarisation. 



base 


NH 3 


H 2 


HF 


STO-3G 


0.703 


0.679 


0.507 


4-31G 


0.905 


1.026 


0.897 


6-31G* 


0.768 


0.876 


0.780 


6-31G** 


0.744 


0.860 


0.776 


limite HF 


0.653 


0.785 


0.716 


Exp. 


0.579 


0.728 


0.716 



Tab. 8.3: Moments dipolaires (u.a.) 



La metnode de coys cernarai (counterpoisej est tres utilisee pour estimer 1 erreur ae superposition. 
L'energie du complexe et des monomeres est calculee avec la meme base AU B, c'est a dire dans la base 
du complexe. 

EfjE AB - E A (A U B) - E B (A U B) 

8.1.7 Pseudopotentiels de cceur 

Tous les electrons d'un atome ne jouent pas le meme role. Ceux des couches internes (electrons de 
cceur) ne participent pas aux liaisons chimiques, alors que les electrons de la couche de valence sont 
les plus actifs. II est avantageux de remplacer les cceurs des atomes par des potentiels effectifs. En ne 
traitant explicitement que les electrons de valence on ne perd pratiquement aucune information sur les 
proprietes physico-chimiques des molecules, mais Ton reduit de facon tres significative le volume des 
calculs a effectuer, surtout si le systeme etudie contient des atomes lourds. 

Dans le formalisme des pseudopotentiels de cceur les electrons des couches internes sont simules 
par un operateur monoelectronique appele pseudopotentiel et qui est extrait en inversant l'equation de 
Schrodinger atomique. 

On distingue trois types de pseudopotentiels de cceur : 

- Les pseudopotentiels locaux, c'est a dire de la forme : 

W PS {r) 

- Les pseudopotentiels semi-locaux : 

„ Imax 

w PS = E w E i y ™)( y ™i 

1=0 m=-l 

- Les pseudopotentiels non locaux : 

W PS = ~^f- + E Mr, ^ 0)><x(r, ^ <f>)\ 

k 

Parmi les pseudopotentiels de cceur les plus utilises je citerai ceux du groupe de Toulouse (Durand et 
Barthelat), ceux de Preuss, Stoll, Igel-Mann, ceux de Hay et Wadt, ceux de Huzinaga et sakai et ceux de 
Bachelet, Schliiter et Hamman. 

Un avantage supplementaire des pseudopotentiels est que les effets relativistes peuvent etre introduits 
dans le pseudopotentiel lui-meme et de ce fait un programme moleculaire non relativiste peut etre utilise 
au calcul de molecules contenant des atomes des quatriemes et cinquiemes lignes. 

8.1.8 Methode Hartree-Fock periodique 

Le nombre d'electrons d'un cristal parfait est infini, celui d'une molecule est fini. Pour traiter les 
cristaux, les surfaces periodiques et les polymeres, on profite de la symetrie de translation pour se ramener 
a un nombre fini de particules, celles de la maille elementaire. 

Dans la methode Hartree-Fock periodique les orbitales cristallines sont representees par des combi- 
naisons lineaires de fonctions de Bloch, elles-memes developpees sur les fonctions atomiques de base de 
chaque maille elementaire : 

*(k;r) =E i « i i(k)^(k;r) 

$(k;r) =E„X?(r)exp(kr) (8.32) 

ou x] ( r ) represente la j eme orbitale de la maille caracterisee par le vecteur de translation n du reseau 
direct. Les coefficients ay sont calcules en resolvant pour chaque vecteur k du reseau reciproque de la 
zone de Brillouin : 

F(k)A(k) = E(k)S(k)C(k) (8.33) 



Hin tait on utilise une grille de points K aont la dimension vane avec la nature du reseau. 

Cette methode permet de calculer la structure electronique de systemes periodiques en 1-2-3 dimen- 
sions. Elle est bien adaptee aux systemes covalents et ioniques, les resultats sur les metaux sont en general 
moins bons. 

Cette methode a ete utilisee avec succes pour calculer des surfaces d'energie en fonction des parametres 
cristallographiques, pour etudier les transitions de phase a l'etat solide, pour calculer des constantes 
d'elasticite et des modules de rigidite, pour tracer des cartes de densite electronique et de potentiel 
electrostatique ainsi que pour calculer des profils de Compton et des facteurs de structure. 

8.1.9 Proprietes des fonctions Hartree-Fock 

Theoreme de Koopmans 

Si un systeme quantique est decrit par le determinant Hartree-Fock d'ordre TV 

D = -^=det{y>iv?2 • ■ .(fk-ifkfk+i ■ ■ -<Pn} (8.34) 

La description approchee la plus simple du meme systeme ionise est un determinant d'ordre N — 1 
construit sur les memes spinorbitales moins une : 

Aon = a^D = j det{yiy 2 • ■ • • -Vjv} (8.35) 

- 1)! 

ou est un operateur d'annihilation. Le potentiel d'ionisation I p est egal a la difference des energies 
electroniques du systeme apres et avant l'ionisation; dans le cadre de cette approximation : 

Ip » (Aon |^| Aon) " (D\H\D) 

N 

« -(A|h 1 |A>-^{(Aj|Aj>-(*J'|j*>} 

3 = 1 

rs -(k\F\k) = -e k (8.36) 

Cette expression approchee des potentiels d'ionisation permet de donner une signification physique aux 
energies orbital aires. 

Theoreme de Brillouin 

On appele determinant monoeexcite du determinant D, le determinant construit sur les memes 
spinorbitales a l'exception de la spinorbitale <fh remplacee par une spinorbitale virtuelle ip^. 

D t = -^=^{fif2---f t i---fN} (8.37) 

ou encore 

Dl = a- k a\D (8.38) 

dans cette equation la substitution de la spinorbitale <fh par la spinorbitale ip^ est exprimee a l'aide des 
operateurs d'annihilation et de creation a+ . 

L'application des regies de Slater permet de calculer l'element de matrice de l'operateur H : 

N 

{D%\H\D) = (n^k) + ^2(jjii\ki) - (/ii|iA) 
= (l*\F\k) 

= (8.39) 



ineoreme d Mellmann- J? eynman 

Le theoreme d'hellmann-Feynman permet de calculer la derivee premiere de l'energie par rapport a 
un parametre a dont depend l'hamiltonien (coordonnees nucleaires, champs externes, etc.), pour une 
fonction d'onde exacte : 

^ = <*l^l*> (8-40) 
da oa 

Cette relation reste valable pour une fonction d'onde Hartree-Fock en effet : 

*w%-m + {¥ww + w\¥) (8.41) 

da da da da 

les deux derniers termes correspondent a une variation infinitesimale de la fonction d'onde Hartree-Fock : 

8{D\H\D) = {SD\H\D) + {D\H\5D) (8.42) 

quantite nulle par definition des spinorbitales Hartree-Fock, done : 

'f = iD\f\D) ( 8 .4 3 ) 
da oa 

Mat rice densite et densite electronique 

La fonction 

N 

t=i 

est appelee matrice densite de Fock-Dirac. Cette fonction est invariante par transformation lineaire des 
spinorbitales. Les spinorbitales etant orthonormees on a : 

i 3 
i 3 

= ^2 l p i (x) l p*(x')=p(x,x') (8.45) 



et 



N 



t=i ' 



N 



J p(x,x)dx = J2 J <fii(x)<p*(x)dx = N (8.46) 
La densite electronique correspond a : 

p{v) = Y J Vi{*)ti{*) (8-47) 

t=i 

Pour un systeme a couches fermees son developpement sur les fonctions de base est donne par : 

occ 

P( r ) = ^E^^^W^^W = ^2 p ^xA T )Xnu(r) (8.48) 



Populations ae iviulliKen 

En developpant l'equation 8.46 sur les fonctions de base il vient : 

J p(r)dr = Y, P*» J X^)Xnu{v)dv = ^(PS)„ = N (8.49) 

la matrice PS est interpretee comme etant la representation de la densite electronique sur la base. Ses 
elements peuvent etre regroupes en fonctions des atonies sur lesquels les fonctions de base sont centrees : 

X)( ps )/'/' = EE( ps )/'/' = E^ = JV ( 8 - 5 °) 

It A pEA A 

Qa est la population electronique de l'atome A. La charge nette de Mulliken est donnee par la relation : 

q A = Z A - Q A (8.51) 

Les charges de Mulliken ne sont definies qu'en fonction de la base utilisee. Leur valeur depend done de la 
facon dont le calcul a ete effectue, ainsi a deux densites identiques, l'une obtenue par un developpement 
monocentrique et l'autre par un developpement multicentrique correspondent des charges de Mulliken 
totalement differentes. Les charges de Mulliken n'ont aucune signification physique intrinseque. 



Chapitre 9 

Les Methodes semiempiriques 



9.1 Methodes autocoherentes 

La resolution des equations Hatree-Fock par la technique de developpement sur une base de fonctions 
connues suppose tout d'abord la construction de la matrice de Fock dont l'expression pour les systemes 
a couches fermees: 



requiert revaluation de nombreuses integrales mono et bielectroniques. En particulier, pour une base 
de dimension N on doit evaluer N(N + l)(iV 2 + N + 2)/8 integrales bielectroniques. D'autre par la 
complexity algorythmique de la construction de la matrice de Fock est proportionelle a N 3 . II est done 
utile de disposer de methodes approchees permettant d'obtenir a moindre cout de calcul des resultats 
semi-quantitatifs. Diverses methodes semi empiriques ont ete developpees dans ce but. Nous citerons les 
methodes CNDO, CNDO/2, INDO, MINDO, MNDO et la methode de Hiickel. 

9.2 Un exemple de methode approchee : la methode de Hiickel 

9.2.1 Description de la methode. 

La methode de Hiickel s'applique aux systemes moleculaires organiques possedant un plan de symetrie 
<r, e'est a dire pour lesquels C s est un sous groupe. Les orbitales, dans ce cas sont soit de symetrie A' soit 
de symetrie A". L'operateur de Fock est totalement symetrique et par suite, les elements de la matrice 
de Fock entre fonctions de symetrie differentes sont nuls. La diagonalisation de la matrice de Fock peut 
alors factorisee par blocs. Dans la methode de Hiickel on ne considere que les fonctions de symetrie A" , 
ce qui correspond au systeme it des molecules. La base est constitue d'une orbitale p z par atome. D'autre 
part on neglige le recouvrement entre orbitales differentes, ce qui correspond a choisir 




(9.1) 



S = 1 



(9.2) 



on doit done resoudre F equation au valeurs propres 



HC = eC 



(9.3) 



Les elements de la matrice H ont une forme extremement simplifiee: 



H. 



a x 



H 



H, 



Pxc 




atomes adjacents 



smon 



(9.4) 
(9.5) 
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rour le car bone 



pour un heteroatome (X=N,0) 



a c 
Pec 



-U.26eV 
-2.5eV 



a x = etc + h x Pcc 
Pxc = x Pec 



(9.6) 
(9.7) 



(9.8) 
(9.9) 





h x 


kcx 


=0 


1.0 


1.0 


-0 


2.0 


0.8 


=N 


0.5 


1.0 


-N 


1.5 


0.8 



Les orbitales molecul aires de symetrie it s'expriment en fonction des orbital es de base (2p z ) Xy.'- 



ifiXfi 



L'energie totale du systeme it est: 



(9.10) 



(9.11) 



oil ri; represente l'occupation de la ieme orbitale et son energie. La charge de l'atome ji (oil est centree 
la fonction de base es ^ donnee par la relation: 



i = l 

enfin l'indice de liaison it entre les atonies indices [i et v est: 

^ nocc 



(9.12) 



(9.13) 



i = l 



La methode de Hiickel est tres interessante pour etudier les systemes aromatiques pour lesquels on definit 
l'energie de delocalisation: 

Edeioc. = E T - m x (2a + 2/3) (9.14) 

to est le nombre de double liaison et 2a + 2/3 est l'energie d'une liaison it dans le cas d'un systeme sans 
delocalisation comme l'ethylene. Dans le cas d'un radical on compte a + /3 par electron ir celibataire. 

9.2.2 Exemple d'application : polyenes. 

Dans le cas des polyenes alternes, la matrice de Hiikel 'ecrit: 



H 



a 


P 





. 


. 





P 


a 


p 


. 


. 








P 


a 


p • 


. 




























. 


• p 


a 



(9.15) 



en laisant le cnangement ae variable x — j 



detH = p n D n 



H 



a; 1 
1 a; 1 
1 a: 1 





Les determinant D„ satisfont la relation de recurence: 











1 X 



D n = xD n _i - D n _ 2 
Pour le butadiene (n=4), les racines de l'equation seculaire sont 

Xl = 1(1 + V5) 

X'2 



X3 
X\ 



les energies des orbitales occupees sont: 



ti 

«2 



= "^(1-V5) 



a +-(1 + Vb)p 



L'energie totale est donnee par 
et l'energie de delocalisation par 



E T = 4a + 2VE/3 
E de i oc = (2^5 - 4)/? = 0.47/3 



Chapitre 10 



Au dela de Papproximation 
Hartree-Fock 

10.1 La correlation electronique 

Nous avons deja montre quelles etaient les limites de la methode Hartree-Fock, nous allons maintenant 
examiner les strategies qui permettent d'ameliorer la fonction d'onde en prenant l'approximation Hartree- 
Fock comme point de depart. 

L'energie de correlation est definie par : 

-Ecorr = -E'exact-NR ~ EjiF 

C'est la quantite importante de toutes les methodes post-SCF. Deux strategies peuvent etre utilisees pour 
aller au dela de l'approximation Hartree-Fock et calculer l'energie de correlation : 

1. Developper la fonction d'onde sur plusieurs determinants : 

I*) = Yl CkDk 

K 

C'est la technique employee par les methodes IC, MCSCF, CASSCF, MP/2, MP/4. 

2. Premultiplier la fonction Hartree-Fock par un prefacteur de Jastrow qui a pour fonction de correler 
les particules et de creer le trou de Coulomb : 

|¥)raexp(^/(ry)) x \V HF ) 

ij 

Cette approche est utilisee dans la methode Monte-Carlo variationnelle (VQMC), dans la fonction- 
nelle de Colle et Salvetti et dans la methode des clusters couples. 

10.2 Developpement sur une base de determinants 

A partir des equations Hartree-Fock on peut construire deux types de spinorbitale : 
les TV spinorbitales occupees qui sont utilisees pour la fonction I^jjf} 

les spinorbitales virtuelles orthogonales aux spinorbitales occupees, il en existe une infinite. 
L'ensemble des spinorbitales occupees et virtuelles forment une base complete parce que l'operateur de 
Fock est hermitique. A partir du determinant Hartree-Fock Do on peut construire des determinants 
monoexcites Df, diexcites D'j, triexcites -D/*Jf, quadriexcites Df-^ etc., dans lesquels une spinorbitale 
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occupee i est remplacee par line spmorDitale virtuelle s. Oes determinants sont symDolises par leur mitiale 

■S D /l . Q . 

Les differentes methodes consistent, avec des techniques differentes, a calculer les coefficients du de- 
veloppement. 

10.2.1 MCSCF et CASSCF 

Dans ces deux methodes on utilise un tres petit nombre de determinants. En contre partie, les spi- 
norbitales sont optimisees simultanement au calcul des coefficients du developpement. 

Dans la methode MCSCF, les determinants sont choisis par l'utilisateur en fonction de criteres phy- 
siques et de son experience. Dans la methode CASSCF (Complete Active Space Self Consistent Field) on 
effectue d'abord une partition des spinorbitales en trois groupes : 

- les spinorbitales de cceur, 

- les spinorbitales actives qui comprennent les spinorbitales de valence ainsi que quelques spinorbitales 
virtuelles jugees importantes, 

- les autres spinorbitales virtuelles. 

On genere ensuite tous les determinants possibles a partir des spinorbitales actives et en gelant les orbitales 
de cceur. Le calcul MCSCF est ensuite effectue sur ces determinants. 



10.2.2 Interaction de configuration variationnelle 



On appelle configuration un groupe de determinants correspondant a une valeur propre degeneree de 
l'operateur de Fock. La fonction d'onde etant representee par le developpement : 



K 

oil ) es t une configuration, le calcul variationnel des coefficients du developpement revient a diago- 
naliser la matrice H dont les elements sont : 

H IJ = {* I \H\*j) 

Un grand nombre de ces elements de matrice sont nuls. Les regies de Slater permettent de determiner si 
un determinant est nul ou non. 

- les elements de matrice entre determinants qui different par plus de deux spinorbitales sont tous 
nuls, 

- les configurations monoexcitees n'interagissent pas avec le determinant de reference, 

- les elements de matrice entre determinants de spins differents sont nuls. 
La matrice H se presente sous la forme : 

($o|#|$o) {$o\H\D) 

{S\H\S) (S\H\D) (S\H\T) 

(D\H\* ) (D)\H\S) (D\H\D) (D\H\T) (D\H\Q) 

(T|^|5> {T\H\D) {T\H\T) \t\H\Q) 

(Q\H\D) (Q\H\T) (Q\H\Q) 



H = 



(10.1) 



Si Ton utilise une base complete de configurations, la valeur propre la plus basse correspond a une borne 
superieure de l'energie non relativiste de l'etat fondamental, les valeurs propres suivantes a des bornes 
superieures des etats excites. 

Les coefficients preponderants du developpement sont d'abord celui de la configuration de reference, 
puis ceux qui correspondent aux diexcitations. 

Dans la pratique, on utilise un nombre fini d'orbitales virtuelles qui est determine par la dimension de la 
base. Beaucoup de calculs sont limites aux seules mono et diexcitations (CISD). Si les configurations sont 



generees a partir de plusieurs conngurations de reierence (ou parentsj, le calcul est appele multireierenciel 
(MRDCI, MRDSDCI). Dans les calculs MRDCI on utilise les diexcitations des configurations de reference, 
dans le MRSDCI les mono et diexcitations. 

10.2.3 Approche perturbative 

Le calcul des coefficients du developpement de la fonction d'onde peut etre alternativement effectue 
a l'aide d'une methode de perturbation. La partition de l'liamiltonien, en un hamiltonien d'ordre zero et 
un potentiel de perturbation, la plus utilisee est celle de M0ller et Plesset : 



L'liamiltonien d'ordre zero est l'liamiltonien Hartree-Fock dont on sait calculer les etats propres et le 
potentiel de perturbation est donne par : 



La fonction d'onde perturbee au premier ordre ne fait intervenir que les configurations diexcitees par 
rapport a la configuration de reference. Les monoexcitations sont prises en compte au troisieme ordre de 
perturbation, les tri et quadriexcitations n'interviennent qu'au quatrieme ordre. Selon l'ordre de pertur- 
bation les calculs sont designes MP/2, MP/3 et MP/4. 

10.3 Fonctions correlees 

La correlation electronique faisant apparaitre le trou de Coulomb peut etre introduite dans les fonc- 
tions de base elles memes. C'est possible pour les systemes a deux electrons mais difficilement generali- 
sable. Pour He les fonctions d'onde de Hylleraas et de Pekeris sont de ce type ainsi que les fonctions de 
James et Coolidge et de Kolos et Wolniewicz pour la molecule d'hydrogene. 

Pour les systemes polyelectroniques il est plus avantageux de multiplier le determinant de Slater (HF) 
par une fonction des distances entre electrons appelee prefacteur de Jastrow : 



Les fonctions de correlation /(r,j) tendent vers — oo quand r,j tend vers zero, d'autre part, afin de ne pas 
compter deux fois le trou de Fermi pour des electrons de spins antiparalleles, ces fonctions sont afiectees 
d'un coefficient suivant la nature des spins. 

Ce formalisme se prete mal a une resolution directe d 'equations integro-differentielles et differentes 
techniques peuvent etre utilisees pour optimiser les parametres du prefacteur de Jastrow. 

10.3.1 Monte-Carlo variationnel (VQMC) 

Dans cette methode, pour chaque configuration {i} des particules, ici configuration signifie ensemble 
de positions instantanees, on definit une grandeur appelee energie locale : 



-^exact — H HF ^ 



(10.2) 




(10.3) 



Y[exp(f( rij )) 



E({i}) = 



(10.4) 



Pour une fonction d'onde exacte, cette grandeur est constante et egale a la valeur propre de l'operateur H . 
Pour une fonction d'onde approchee, l'energie locale fluctue autour d'une valeur moyenne. La meilleure 



tonction approcliee est celle qui minimise les nuctuations ae 1 energie locale et on est amene a mmimiser 
la fonctionnelle : 



2 _ 
^opt - 



E,;T(^({i})-^({i}) 



,N. 



(10.5) 



qui correspond a la variance de l'energie locale. E est la moyenne de l'energie locale sur l'ensemble des 
iV opt configurations retenues et 



*({.}) 



*o({i}) 



une fonction de ponderation ou \Po({i}) es ^ l a fonction de depart avant minimisation. Les 7V opt configu- 
rations sont engendrees par 1'algorithme de Metropolis. 



10.3.2 Fonctionelle de Colle et Salvetti 

Dans cette methode, le prefacteur de Jastrow ne contient pas de parametres variationnels et l'energie 
de correlation peut etre calculee directement en utilisant une fonctionnelle de la densite Hartree-Fock. 
Ce type de correction est applicable aux systemes periodiques pour lesquels le developpement en serie de 
determinants est impossible. 



10.3.3 Methode des cluster couples 

Dans cette methode qui peut etre egalement presentee comme un developpement en serie de determi- 
nants, le prefacteur de Jastrow a la forme d'un operateur : 

I*} =exp(f)|$) (10.6) 

ou |$) represente la fonction Hartree-Fock et T est un operateur faisant intervenir des operateurs de 
creation et anihilation dans un formalisme trou-particule. Les element de matrices sont evalues a l'aide 
de techniques diagrammatiques. En fait la fonction d'onde peut etre reecrite sous la forme d'un develop- 
pement : 

JV 

|¥)=X)6l*> (10.7) 

t=i 

ou les operateurs C\ qui ont pour expression : 

Ci = fi 
C 2 = f 2 + Ml 

c 3 = n + f^ + i^ (io.8) 

generent des configurations mono, di, triexcitees etc. La methode des clusters couples appartient a la fois 
aux techniques de prefacteur et de developpement. 



10.4 Methode de la fonctionnelle de la densite 
10.4.1 Theorie 

En 1964, Hohenberg et Kohn ont montre qu'il existe une fonctionnelle universelle permettant de 
calculer l'energie d'un systeme quantique a N particules a partir de la densite : 



p(r) =*'(r)*(r) 



pour im potentiel donne v[r) 1 energie de 1 etat londamental est : 

E[p\ = J v(r)p(r)dr + F[p\ (10.9) 

Pour un systeme constitue d'electrons et de noyaux, v(r) correspond a l'attraction coulombienne electron- 
noyaux et F[p] a la somme des contributions des energies cinetique, coulombienne electron-electron, 
d'echange et de correlation. II est possible de reecrire cette equation sous la forme : 

E[p] = J v(r)p(r)dr + ± J J + G[p] (10.10) 

avec G[p\=T[p\ + E xc [p\ 

T[p] est la contribution cinetique et £" X c la contribution d'echange-correlation. Pour une densite 
uniforme on montre que : 

^xc[p] = I p{r)e xc (p(r))dr (10.11) 

oil exc(/?( r )) est l'energie d'echange et de correlation par electron appartenant a un gaz d'electron de den- 
site uniforme p. Si la densite presente des fluctuations on developpe £"xc[/3] en serie ce qui fait apparaitre 
des corrections dependant des gradients successifs de la densite. 

Pour un systeme a N electrons, Kohn et Sham ont derive une expression variationnelle qui permet 
de calculer la densite. Si Ton represente la fonction d'onde du systeme par un determinant construit sur 
des spinorbitales orthonormees <f>i, alors : 

N 

p(r) = ^*(r)^(r) (10.12) 

i = l 

L'application du principe variationnel a l'equation 10.10 conduit a un systeme de N equations couplees : 

{_Iv 2 + ^(r) + /JxcMr))}^} = eMi) (10.13) 
similaires aux equations Hartree-Fock et appelees equations de Kohn-Sham, dans cette expression 

<f{r) = v(r) + J j^^dr' (10.14) 

et 

= d^r) ip{r)excip{r))) (1 °- 15) 

Le probleme crucial a resoudre est la determination de la fonctionnelle /Jxc(p(r)) pour laquelle il 
n'existe pas d 'expression analytique. 



10.4.2 Approximation de la densite locale (LDA) 

Bien qu'il n'existe pas de forme analytique de la fonctionnelle ^xc(p), on peut deriver des expressions 
approchees dans le cas du gaz d'electrons de densite uniforme. Cette technique permet de parametrer 
une expression de la fonctionnelle a partir de differentes valeurs de la densite. Dans le cas des sytemes 
reels ou la densite n'est pas uniforme, cette approximation n'est que localement valable. Les corrections 
de gradient permettent de representer correctement les regions de grande fluctuation. 

Dans le cas des systemes de fermions l'approximation LSD (Local Spin Density) permet d'introduire 
la densite de spin en partitionnant la densite totale en une contribution des electrons de spin a et une 
contribution des electrons de spin /?. Les fonctionnelles les plus precises ont ete derivees par Perdew et 
Zunger a partir des calculs Monte-Carlo quantique de Cepperley et Alder sur le gaz d'electrons. Plus 
recemment Perdew a propose une approximation de gradient generalise tres satisfaisante dans les regions 
de grande variation de la densite. 



iu.i.o ivieinoae ae v^ar ex, .rarrineno 



Les methodes de fonctionnelle de la densite sont encore tres peu exploiters par les quantochimistes 
pour effectuer des calculs moleculaires, par contre, en physique, du solide cette methode est beaucoup 
utilisee. Parmi ces methodes, celle de Car et Parrinello succite un grand interet et semble promise a 
beaucoup d'avenir. 

Les orbitales cristallines de Kohn-Sham sont developpees sur une base d'onde planes. La recherche 
des coefficients des fonctions de base n'est pas effectuee par la diagonalisation successive de la matrice de 
l'hamiltonien effectif, mais par une dynamique moleculaire oil ces parametres sont assimiles aux positions 
de pseudo particules de masse Active. Afin de trouver le minimum absolu, cette dynamique est effectuee 
a differentes temperatures fictives afin de simuler un recuit (methode du recuit simule). L'avantage de la 
technique de Car et Parrinello est qu'il est possible d'optimiser la structure, c'est a dire les positions des 
noyaux, simultanement. Pour cela, une dynamique moleculaire est effectuee sur les noyaux dans le champs 
des electrons. Le volume des calculs reste raisonnable parce que les fonctions de base sont independantes 
des positions des noyaux. 
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Annexe A 

Systemes de coordonnees 



A.l Coordonnees orthogonales generalisees 

Les coordonnees orthogonales generalisees {xi, x 2 , £3} sont definies par les vecteurs unitaires ei, e 2 , e3. 
L'element differentiel ds s'exprime, en fonction de dx\,dx 2 ,dx 3 et des coefficients h\,h 2 ,h 3 



ds = eihidxi + e 2 h 2 dx 2 + e 3 h 3 dx 3 



et 



ds 2 = h\d 2 xi + h\d 2 x 2 + h\d 2 x 3 
les coefficients h s sont determines par la relation suivante : 



h s ei = 



dx . dy , dz 
1— +j-T^ + k : 



da 



'da 



da 



done 



/iiei 
h 2 e 2 
/iiei 



. dx dy 

h Ja hk ^ — 

dxi ox\ ox\ 



dz 
dx 
dz 



. dx . dy , 

dx 2 ^ dx 2 dx 2 

. dx .9 dz 
1 o VS-dx?, + k- — 

dx 3 y dx 3 



comme ei,e2 et e3 sont orthogonaux 



v da 



dz 
'dx. 



(A.l) 
(A.2) 

(A.3) 

(A.4) 
(A.5) 
(A.6) 

(A.7) 



A. 1.1 Expression des operateurs vectoriels. 
gradient 



ei du e 2 du e 3 du 
hi dxi h 2 dx 2 /13 dx$ 



divergence 



(A.8) 



V • V = 



1 



hih 2 h 3 



(A.9) 
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laplacien 



V 2 u 



1 



hih 2 h 3 



d h 2 h 3 d d h x h 3 d d h x h 2 d ' 

dxi hi dxi dx 2 h 2 dx 2 dx 3 h 3 dx 3 



(A.10) 



rotationnel 



V A V = 



1 



hih 2 h 3 



/iiei h 2 e 2 h 3 e 3 
d d d 

dx\ dx 2 dx 3 

hiVi h 2 V 2 h 3 V 3 



element d'integration 



(A.11) 



dr — h\h 2 h 3 dx\dx 2 dx 3 
A. 1.2 Coordonnees cylindriques 



(A.12) 



x = r cos 



y — r sin ( 



(A.13) 



z^z 
2 



hee e 

h z e 

h . 

hi 

hi 
h r 

he 
K 



V« = - 



V 2 u 



i cos + j sin 6* 
— ir sin + jr cos 6* 
k 

cos 2 + sin 2 = 1 
r 2 sin 2 + r 2 cos 2 = 
1 
1 
r 
1 
1 



du 1 3w 3w 

r dr r 89 z dz 



15 3 1 d 2 d 2 
r dr dr r 2 dO 2 dz 2 



dr = rdrdOdz 



(A.14) 
(A.15) 
(A.16) 
(A.17) 
(A.18) 
(A. 19) 
(A.20) 
(A.21) 
(A.22) 

(A.23) 

(A.24) 
(A.25) 
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X 




v sin cos <b 




(A. 26) 
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Annexe B 

Integrates utiles 



In{p) — J x n exp (— px)dx = n\p~ n ~ l 



A n (p) — / x n exp (-px)dx = nip'"' 1 exp (-p) ^ — 

{ k = 

1 

B n (rp) = j x n exp(-T P x)dx = -A n {rp) - (-l) n A n (-T P ) 



(B.l) 
(B.2) 
(B.3) 



105 



Index 



S 2 , 55 
S 2 , 55 
S, 55 

elements, 9 

equation de Schrodinger, 19 

Abbeg, R., 14 
Anaximandre de Milet, 10 
Aristote, 10 

Bergman, T., 11 
Berthollet, C, 12 
Berzelius, J. J., 12, 13 
Bohr, N., 14 
Boscovitch, R., 11 
boson, 21 
Buffon, G., 11 
Burrau, 0., 15 

Colle, R., 90 
Condon, E. U., 15 
constante du mouvement, 22 
Coolidge, A. S., 15, 69 
Couper, A. S., 13 

Democrite, 10 
Dalton, J., 12 
Daudel, R., 16 
Davy, H., 13 
de Broglie, L., 15 
Descartes, R., 11 
dihydrogene, 63 

effet Stark, 26 
Ehrenfest, 21 
Empedocle, 10 

Faraday, M., 13 
fermion, 21 

fonction de Guillemin et Zener, 49 
fonction de James, 50 
fonction de James et Coolidge, 69 
fonction de spin, 53 
fonctionnelle de la densite, 90 



Fontenelle, J.-B., 9 
Fourcroy, F., 12 
Frankland, E., 13 

Gay-Lussac, J.-L., 12 
Geoffroy, C.-J., 11 
Guillemin, V., 49 
Guyton de Morveau, L. B., 12 

H+,45 

hamiltonien, 19 
Heisenberg, 21 
Heitler, W., 15, 63 
Helmholtz, H., 13 
Helmont, J. B. van, 11 
Hohenberg, 90 
Hylleraas, 15 

indiscernabilite, 20 
interaction de configuration, 68 

James, H. M., 15, 50, 69 

Kant, E., 12 
Kekule, A., 13 
Kohn, W., 90 
Kossel, W., 14 

Lavoisier, A. L., 12 
Le Bel, J. A., 13 
Lewis, G. N., 13, 14 
liaison de valence, 68 
London, F., 15, 63 
Loschmidt, J., 13 

mecanique quantique, 15 

methode de variation des constantes, 26 

methode de variation perturbation, 27 

Macquer, J., 12 

Mendeleev, D., 13 

Millikan, R. A., 13 

Newton, I., 11 

operateur energie cinetique, 19 



106 



operateur energie potentielle, iy 
operateur de spin, 53 
orbitale mpleculaire, 68 

Paracelse, P. A., 10 
Pasteur, L., 13 
Pauli, W., 15, 53 
Pauling, L., 10, 15 
permutation, 21 

perturbation dependant du temps, 26 
perturbation statique, 23 
Platon, 10 

prefacteur de Jastrow, 90 
principe d'antisymetrie, 21 
principe de Pauli, 15, 53 
Prout, W., 12 

regie de l'octet, 14 
Ramsay, W., 14 
relation d'Heisenberg, 21 
relations d'Ehrenfest, 21 
Rutherford, E., 14 
Rydberg, J. R., 12 

Salvetti, O., 90 
Schrodinger, E., 15 
Stahl, G. E., 12 
Stark, J., 14 

theoreme d'Hellmann-Feynman, 23 
theoreme du viriel, 22 
theorie de la valence, 14 
theorie des perturbations, 23 
theorie dualistique, 13 
Thales de Milet, 10 
Thomson, J. J., 13 
trou de Fermi, 54 

Van't Hoff, J. H., 13 

Werner, A., 13 
Wilson, H. A., 13 



Zener, C, 49 



